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1 Introduction 

The classical integrable Neumann system on an n-dimensional sphere is known to 
be a toy model for the Jacobi-Mumford systems introduced in [42] and describing 
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translationally invariant flows on the Jacobians of hyperelliptic curves and repre¬ 
sented by 2 X 2 Lax pairs. 

During the last two decades various generalizations of the JacobDMumford 
systems associated with families of r-gonal and generic curves were constructed 
([isiEaiMiEni). The paper [2S] also introduced such systems linearized on special 
types of Prym varieties isomorphic to Jacobians. All such generalizations are closely 
related to integrable flows on coadjoint orbits of the loop algebra gl(r), with Lax 
matrices being rank r perturbations of a constant matrix. The complete algebraic 
geometric description of such flows in the generic case was given in the series of 
papers [HI [Dill [3]. 

There have also been several important results on discretizations (Backlund 
transformations) of the above systems, which are described by translations on the 
Jacobians of hyperelliptic, r-gonal, and generic curves, and which can be related to 
addition theorems for meromorphic functions on such Jacobians ( |561136t I21|). 

The present paper contributes to these fields of interest in two ways. 

First, as a further generalization of the Jacobi-Mumford systems, we consider 
natural analogs of the Neumann system on the Stiefel variety Vn,r < n), the set 
of n X r matrices X satisfying the constraints X'^X = 1^. These analogs were first 
introduced by Reiman and Semenov m by presenting a ’’big” (n x n) Lax pair, 
which is closely related to the Lax pair for the symmetric Clebsch-Perelomov rigid 
body problem [35]. For r = 1 it had been given by Moser in m- It also appears 
that for r > 1 the big Lax pair itself does not define yet the Neumann system on Vn,r 
uniquely. Namely, there exists a family of integrable Neumann flows with different 
SO(n)-invariant metrics on Vn,r that possess the same big Lax representation. 

Using the previous results of [33|, in Section 3 we show that these analogs of 
the Neumann systems on Vn^r also admit ’’small” (2r x 2r) Lax representations, all 
having the symplectic block form 

|L(A) = [L(A),ArjA)], 

^^(Aln - X^(AI„ - A 

[lr-P^{Xln-A)-^P -P^{Xln-A)-^Xj ' 

where A = diag(ai,..., a,i), P is the n x r-momentum satisfying the constraint 
X^P + P"^X = 0, A is a rational spectral parameter, and the matrix which is 
not shown here, depends on a parameter k of the metric. 

The matrix L(A) can be regarded as a direct generalization of the 2x2 Lax 
matrix for the classical Neumann system f |42l [5l [30]i and fits to the already studied 
class of rank d perturbations of the constant matrix A. The latter are described by 
the Lax matrices of the form 


n . r 

+ Egl(d), 

where Y is constant and A/) depend on the variables of the corresponding dynamical 
system. However our L(A) belongs to the symplectic loop subalgebra s'p(2r) C gl(2r) 
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and has a quite specific structure, moreover, V is not diagonalizable, it contains 
a Jordan block. This causes the spectral curve S of L(\) to have extra strong 
singularities at the infinity. So, the previous results of mw concerning the genus 
of S, number of its infinite points, etc, need to be modified in order to treat this 
case. 

We show that the Neumann systems on are integrable in a non-commutative 
sense m), and calculate the dimension 5 of its generic minimal invariant tori. 

Calculating the order of singularity and the genus of the regularized curve S, 
taking into account the involution a : S —>■ S, we then demonstrate that, up to the 
action of a discrete group of reflections, generic complex invariant manifolds of the 
Neumann systems on T*Vn,r are open subsets of affine Prym varieties Prym(5, a) of 
the same dimension 6. The latter are algebraic noncompact subgroups of generalized 
Jacobian varieties Jac(5,oo), the extensions of Jac(5) by 

It should be noted that affine Prym varieties appear as complex invariant mani¬ 
folds in several classical integrable systems on Lie groups before their reduction, for 
example, the spatial rigid body motion in the Clebsch case of the Kirchoff equations 
(see e.g., m)- However, to our knowledge, in the existing literature such manifolds 
were not given the corresponding algebraic description. 

Next, the Marsden-Weinstein reduction of the Neumann systems on T*Vn,r by 
the group SO{r) are also integrable, and their complex invariant tori are shown to 
be usual Prym varieties Prym(5,(T) C Jac(5) (again, up to the action of a discrete 
group). The same holds for the reductions of the systems on the cotangent bnndle 
of the Grassmann variety Gn,r: which are integrable in the commutative sense. 

Second, we construct a family of (multi-valued) integrable discretizations (Backlund 
transformations) IB,. : {X,P) —>■ {X,P) of the Neumann system on T*Vn,r- The 
family is parameterized by A* G C and described by the 2r x 2r intertwining relation 
(discrete Lax pair) 

Z(A)M(A.A.) = M(A.A.)L(A), M(A. A.) = 

where L{\) depends on X,P in the same way as L(A) depends on X,P, whereas 
r(A*) depends on X,X in a symmetric way. Again, the above relation is a direct 
matrix generalization of the 2x2 discrete Lax pair describing a Backlund transfor¬ 
mation of the classical Neumann system (see |29l l36j L The Lagrangian description 
of the map is closely related to a discretization of the Neumann systems on 14,,r 
introdnced in [51] in connection with factorization of matrix operators. 

Each discretization map 53^ preserves the first integrals of the continnons system 
and geometrically is described by translations T on Prym(5,(T). The translations 
are written explicitly in terms of r non-involutive points on the spectral curve S over 
the coordinate A* (i.e., eigenvalues of L(A*)). A choice of such points (partition) 
fixes the branch of the map and the corresponding translation T. Iterations of 
all possible translations generate a lattice of rank on the universal covering 
of Prym(5,iT). The above construction is similar to that used in |41] to describe 
the integrable discrete Frahm-Manakov top on SO{n) and is actually based on the 
procedure of solving a matrix quadratic equation presented earlier in m- 
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A similar matrix generalization can be made for the another discretization of 
the Neumann system on proposed in [39], and which can be considered as an 
example of a generalization of the even Jacobi-Mumford systems. 

In Conclusion we discuss further possible generalizations. 


2 The Jacobi—Mumford and the Neumann systems 

2.1 The standard (odd order) Jacobi—Mumford systems 

Let r = T{R) be a smooth hyperelliptic curve of genus g, whose affine part is given 
by the equation ijp = R{X), where R{X) is a monic polynomial of degree 2g + 1. We 
regard L as a compact Riemann surface having one infinite point oo. Consider a 
generic divisor of g points Pi = ■ ■ ■ ,Pg = {Xg,fig) on L. Following Mum¬ 

ford |42] . the curve and the divisor on it can be associated to three polynomials 
U (A), V (A), and W (A), A G C such that 


U{Xi) = 0, V{Xi)=Hi, i = l,...,g, (2.1) 

W{X)U{X)+ V\X) = R{X). (2.2) 


Hence, the degrees of U (A) and V (A) are at least g and g — 1 respectively. For a 
fixed curve F given by P(A), the polynomial 1F(A) is uniquely determined from the 
condition Then, if U{X) is monic of degree g, then W{X) is monic and has 

degree <7 + 1. That is, 


U (A) — (A — Ai) • • • (A — Xg') — A^ + u\X^ + • • • + Ug, 
= '^l^kY[ V—^ = viX<^~^ + 


k=l j^k 


Xk Xj 


+ Vg, 


(2.3) 


VF(A) = A®+^ + wqX^ -I- wiX^ ^ H-h % = (A - i^i) • • • (A - I'g+i). 


Now let £g be a variety of all non-constant coefficients of P(A), R(A), 1F(A). 
As shown in [32], for a fixed P(A), formula (12.2p defines a set of equations on the 
coefficients which provides a purely algebraic description of an affine part of Jac(F), 
the Jacobian variety of the curve F. The coefficients themselves are meromorphic 
functions on the Jacobian. 

The variety £g itself can be completed to the fiber bundle £g —)• TZ over the 
2g + 1-dimensional space TZ spanned by the coefficients of P(A) and parameterizing 
odd order hyperelliptic curves F = r(P) of genus g (the base), with fibers of £g 
being the Jacobians of the curves. 

Any translationally invariant vector ffow on Jac(r) can be extended (in different 
ways) to a vector ffow on the whole space £g, which leaves the fibers invariant. 
The latter flow is described as an integrable system of differential equations on the 
coefficients of the three polynomials, which is referred to as a Jacobi-Mumford or 
an (odd order) master system (see e.g., [32l [55l [26] i. 


4 



Let i? G 7^, r = r(i?) be the associated hyperelliptic curve and let ^ : L — 
Jac(r) be the Abel mapping given by the following integral with the basepoint oo 


fP 

A{P) = / {uJl, . . . ,UJg)'^, PgL, 

J OO 


where wi,..., Wg are independent holomorphic differentials on the curve F. 

Consider a translationally invariant flow on Jac(r) which is tangent to the image 
of the curve ^(F) C Jac(F) at a finite point A{X*, fi*). Then, fixing A* G C, but not 
R &TI, we get a vector flow on Eg or on Eg. 

The flow can be represented in the following Lax form with A as a rational 
parameter (see my- 


L(A) 

L(A) 

N{X) 


[L(A),iV(A)] 

(V{X) U{X) \ 

Vf^(A) -V{X))^ 

1 f -V{X*) -U{X*)\ 

2(A - A*) V-W^(A*) - (A - X*)U{X*) V{X*) ) 

L(A*) 1/^0 OA 

2(X-X*) 2 [-U(X*) oj' 


(2.4) 


Concerning the Lax equation (j2.4jl . the coefficients of the spectral curve 
det(L(A) - fih) V^(X) - W(X)U(X) = - R(X) = 0 


are the first integrals of the system, and we recover the family of hyperelliptic curves 
F = r(R). When A* tends to oo, after an appropriate time rescaling, the limit flow 
is described by (12.4p with 


=(a 

Remark 1 To a fixed curve F and the given divisor Pi + • • • + on it, one can 
associate the dual divisor Pi + • • • + Rg+i given by the zeros of the meromorphic 
function /(A,/x) = (/r + V{X))/U{X). In view of (j2.1[) and (j2.3[l . f{X,fj,) have simple 
poles at the original divisor and oo, which implies 

-T(Pi) + • • • + A{Pg) + ^(oo) = ^(Pi) + • • • + A{Rg) + A{Rg-\-i). (2-6) 

As follows from (12.211 . one can also write / = W{X)/{fi — I/(A)). Hence, the coordi¬ 
nates Uk = A(Pfc), k = 1,..., g + 1 are the roots of W{X). 

The Neumann system on S^~^. The best known example of an integrable prob¬ 
lem associated to the Jacobi-Mumford system is the Neumann system describing 
the motion of a point on the unit sphere = {{q, g) = 1} C M"', q = {qi ,..., qn) 
with the quadratic potential V{q) = —^{q,Aq), A = diag(ai,... ,an) (see [441 [30] L 
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Upon introducing momentum p = {pi,... ,pn) = q and imposing the constraint 
{p, q) = 0, the motion is described by the equations 


q = p, p = Aq + uq, u =-{p,p) - {q, Aq), i = 
which possesses the 2x2 Lax representation 


E n qiPi spn Qj 

2=1 X—ai Z-/2=l X—ai 

1 _ pj _ qiPi 

Z-/2=l X—ai Z^i=l X—ai, 
0 i\ 

\ + v{p,q) oy 


L(A) = [L(A),iV(A)], L{\) = 

N{\) = 

Let a(A) = (A — ai) • • • (A — a„). For the polynomial Lax matrix 

L(A) = a{X)L[X) 

the characteristic equation |L(A) — pl 2 \ = 0 has the form 

iiiPj - QjPif 




E 




. . (A ci^)(A dj) . A Qi 

and gives the odd order hyperelliptic curve F of genus g = n—\ 

= R{\) = o(A) (A - Cl) • • • (A - Cn-i). 


(2.7) 

( 2 . 8 ) 
(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


The zeros ci,..., c^-i of i?(A) are real (see Lemma 4.5 in Chapter III of [32]) and 
they represent commuting integrals of the Neumann system. 

One can then identify the 2x2 polynomial matrices (I2.10p and (12.41) by setting 
g = n — 1, which gives 

n 2 

U{X) = a{X) ^ = AS + (- tr ^ + {q, Aq))X3-^ + • • • , 

i=i 

n 

V (A) = o(A) = {p, Aq)X^~^ + ■■■ , 

n 2 \ 

1-=AS+i-(tr^ + (p,p))AS + ... , 
which, in particular, implies 

ui =-ti A + {q,Aq), vi = {p,Aq), wq =-ti A-{p,p). (2.13) 


VF(A) = a(A) 


Then the second matrix (1^ coincides with ()2.5p . Hence the Neumann flow is 
linearized on the Jacobian of T, on which it is tangent to Al(r) C Jac(r) at its 
infinite point. 
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The above relations also give the parametrization of in terms of Ai,..., Ag, 
the A-coordinates of the points Pi,...,Pg on T, which now play the role of the 
spheroconical coordinates on 5'"'“^: 


<li = 


{Qji Ai) • • • (cij An—l) 

~ %') 


i = 1,..., n 


(2.14) 


In addition, pj 


can be expressed in terms of A-coordinates of the dual divisor 



{aj - t/j) • • • {aj - Vn) 


Since pj, qf^PiQi are linear functions of the coefficients of U,V,W, they are 
meromorphic function on Jac(r). However, the coordinates Pi,qi themselves do not 
have this property. According to [l2], the following theorem holds 


Theorem 2.1 1) Complex invariant manifolds of the Neumann system with the 

constants of motion Ci factorized by the action of the discrete group hlf gener¬ 
ated by reflections {qi,Pi) —>■ {—qi, —Pi), i = 1, ■ ■ ■ ,n are open subsets of the 
Jacobian Jac(r) of the spectral curve 

2) The complex invariant manifolds themselves are open subsets of unramified cov¬ 
erings o/Jac(r) obtained by doubling some of the period vectors o/Jac(r). 


2.2 Discretizations of the Neumann system on T*5'” ^ 

The first integrable discretization of the Neumann system (12.7j) was found in |56] 
by using the approach of Lagrange correspondences and the idea of factorization of 
Lax operators (see m for the details). 

Combining the results of |36l 154] . below we consider the discretization map 

55 : iq,p) ^ iq,p), 

written in the form containing an extra parameter A* G C: 

q = A-^/‘^{X^){'yq+p), p =-A^/^{X^)q + + jp), (2.15) 

wher43 A(A*) = A*In — A = diag(A* — oi,..., A* — a^). 

Multiplying the second vector equation by q and using the first one, as well as 
constraints {q,q) = 1, {q,p) = 0 , one hnds 

7 = (g,A^/ 2 (A*)g). (2.16) 

The equations (j2.15p . (j2.16jl of the map were obtained in [29] in the context 
of Backund transformations. For any A*, it preserves all the first integrals of the 
continuous Neumann system. Since the variable 7 depends on both sets of variables 

^The above notation slightly differs from that of Chap. 21.6 of [^, where one should replace 
and SI by 1/A* and —A, respectively. 
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in a symmetric way, relations ()2.15p . (12.161) describe the map in an implicit form, 
as most of Backund transformations. The dynamics can be also written in terms of 
discrete Euler-Lagrange equations on (e.g., see [M] i 

A^/‘^{X^){q + q) = I3q, 


where q is the position on ^ previous to q and /3 is the multiplier ensuring that 
q, q remain on 




(2.17) 


This form of equations with = J represents the Heisenberg model of a chain 

with classical spins and was originally studied by Veselov [MlEZj. The Heisenberg 
model has a natural formulation in the pseudo-Euclidean spaces as well (see [32]). 

As was shown in [29], up to the action of the group of reflections {pi,qi) —)• 
{—Pi,—qi), the above map is equivalent to the following intertwining relation (dis¬ 
crete Lax pair): 

L(A)M(A|A*) = M(A|A*)L(A), (2.18) 


where L{X),L{X) have the same structure as in ()2.8I) and depend on the “old” 
variables {p,q) and the “new” ones {q,p) respectively, and 


M(A|A,) = 


-7 1 

A - A* -I- 7 ^ -7 


(2.19) 


Indeed, setting in (I2.18p subsequently A = ai,...,an and calculating the matrix 
M(A|A*) L(A) M“^(A|A*)|A=ai) obtains 


% = ± 


IQi + Pi 
A^ CLi 


Pi = A: 


( 7 ^ + ai- X^)qi + 'jPi 
\/ A^ 0>i 




( 2 . 20 ) 


which is equivalent to (|2.15l) if we fix sign above. 

The relations (I2.15P can also be written in the equivalent alternative form 


p = A^/2(A^)g - 7g, p =-A^/'^{X^)q + 


'yq. 


( 2 . 21 ) 


To find the factor 7 as a function of q,p, one applies the condition {q,q) = 1 to 
the first equation (|2.15p and obtains the quadratic equation 


{q,A \X^)q)+ 2{p, A ^{X^)q)-f + {p, A ^(A*)p) - 1 = 0, (2.22) 


which gives 


7 = 


-{p,A ^{X^)q) ± y/{p,A i(A*)g)2 - (g, A ^{X^)q){{p,A i(A*)p)-l) 


{q,A ^{X^)q) 

or, in view of ( 12 . 111 ) . ( 12 . 121 ) . 

— {PiA ^(A*)g)-|-/i(A*)/a(A*) — V(A*)-|-//(A*) 


7 = 


{q,A i(A*)g) 


U{X, 


(2.23) 



























/i(A*) being the coordinate of one of the two points on the curve () 2 . 12 p over A = A*. 

As a result, the map IB given by ()2.15p . (12.161) preserves the same first integrals as 
the continuous Neumann system, and, in view of Theorem l2.ll its complex invariant 
varieties are open subsets of unramified coverings of the Jacobian of the curve T given 
by (| 2 . 12 p . As also follows from (|2.23ll . the map IB, as well as its inverse *B“^, is 
generally 2 -valued. 

Remark 2 Obviously, equations (j2.15p describe a real map when /i(A*) is real, i.e., 
R{\*) > 0. In particular, i?(A*) > 0 when A* > ai,... , 0 ^, 01 ,... ,Cn-i (see (| 2 . 12 p ). 

Remark 3 Using the above formulas, one can prove that the composition of the 
two branches of IB, corresponding to different choices of sign of /r(A*), gives the map 
In Section [ 6 ] an analog of this property for the discrete system 
on Vn,r will be considered. 

Algebraic geometrical interpretation of the map (12.151) . (I2.16P was given in [56] 
(see also |41] ) and, in the context of discretization of the Jacobi-Mumford systems, 
in |36] . Namely, let IB be the reduction of IB under the action of the reflection group 
Z 2 . Hence, following Theorem 12.11 its generic invariant varieties are open subsets of 
Jac(r). Then, according to the above remark, the two branches of IB just represent 
a ’’shift back” and a ’’shift forward” maps. 

Let, as above, A : T —> Jac(r) be the Abel map. In [56] and [36] the following 
theorem was proven. 

Proposition 2.2 The restriction of IB onto Jac(r) is described by translation by 
the vector 

T = A{P)-A{oo)=A{P), P=(A*,;u*). (2.24) 

Notice that, although the translation vector does not depend explicitly on the 
constants of motion, it depends on them via the moduli of the curve T. The above 
formula allows to write explicit solution for the trajectory {q,p)k in terms of hy- 
perelliptic theta-functions associated to T, whose arguments depend linearly on the 
discrete time A: G Z. As follows from Proposition 12.21 the two branches of IB cor¬ 
responding to P = (A*,/r*), (A*, —p*) are just shifts by opposite vectors ±T. This 
observation is consistent with Remark [3| above. 

Special cases A* = Oj. According to formula ()2.24l) . when (A*,/x*) = {aj,0) the 
translation T is just a half-period in Jac(r) and, in view of (|2.23p . the map IB is 
single-valued. As expected, double iteration IB^ gives the same point in the Jacobian, 
on which pf^q^^Piqi are meromorphic. On the other hand, {q,p) = fB^(( 7 ,p) = 

{-q,-p)- 

This can be checked directly: let us set A* = aj in (I2.20p . Since the denominator 
in j-th equations vanishes, for qj to be finite, one has to set 7 = —pjlqj- Then 
the second equation in (I2.20h gives pj = 'yqj, i.e., Pj/qj = 7 , which, under the next 
iteration of IB, implies 7 = — 7 . Using this relation and iterating (12.151) directly we 
get q = —q, p = —p for any p, q. 

In this sense, the original map IB with A* = aj represents the ’’imaginary unit 
map”. 


9 






























Continuous Limit. As also follows from (12.241) . the vector T tends to zero when 
A* —)• oo, which must give the continuous limit of the map (I2.15p . (12.1611 . We shall 
describe this limit in details in Sectional jointly with the discrete Neumann system 
on the Stiefel variety Vn^r- 


Another discretization of the Neumann system and even order Jacobi— 
Mumford systems. In [49] Ragnisco found another integrable discretization of 
the Neumann problem, which has different first integrals in comparison with the 
classical case. It was assumed that the canonical variables {q,p) are subject to the 
constraints {q,q) = 1, {p,q) = 1/2, and the corresponding discretization mapping is 


q= -{Aq-uq-p), p = aq, 
a 


(2.25) 


where u = {q,Aq) — 1, = {Aq — uq — p)"^. If we replace u by the expression 

{q,Aq) — 2{p,q), the mapping can be extended to the mapping on x M”' = 

{(QtP) I {QjQ) = !}• this case the product {p,q) becomes a first integral. 

The extended mapping (12.251) can be written in discrete Lax form 


L^(A)Mr(A) = Mr{X)Lr{X), (2.26) 


where 



1 ’liPi 

2 A—a. 


\—ai 


Qj \ 

\—ai j 

_1 _ QiPi j ’ 

2 Z-yi=l \—ai / 


0 -1 

X- {q,Aq) +2{p,q) 


and Lji{X) depends on {q,p) in the same way as Lr{X) does on {q,p). For the case 
{PjQ) = "^I‘2‘ the discrete Lax representation was obtained in [50] . 

The characteristic equation |a(A)Lij(A) — /r/| = 0 gives us an even order hyper- 
elliptic spectral curve of genus g — 1, 




= <-(T E 


, i<j 


{qiPj - qjPif 

(A - ai){X - aj) 


E 

2 = 1 


PiQj 
X CLi 


--a(A)(A-ci)---(A-Cn), 

(2.27) 


Cl,... c„ being constants of motion. 

Note that, in contrast to the discrete Neumann system given by ()2.18p . (|2.19l) . 
the mapping (|2.25p does not include a parameter, but has an explicit form. In 
fact, as shown in [21], (12.251) can be included to a family of generally 2-valued maps 
parameterized by A* € P. Each map of the family preserves the same integrals and 
describes a shift on the Jacobian of the curve (I2.27P that depends on A* in a manner 
similar to that of Proposition 12.21 These maps can be regarded as discrete models 
of even master systems considered in [55] [36| . 
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3 The Neumann systems on the Stiefel variety Vn^r 

Recall that the Stiefel variety is the variety of ordered sets of r orthogonal unit 
vectors ei,..., e,. in the Euclidean space (M”', (•, •)), or, equivalently, the set of n x r 
matrices X = (ei • • • e^) satisfying the constraints X = I,.. Thus Vn,r is a smooth 
subvariety of dimension rn — r(r +1)/2 in the space of n x r real matrices (1^) = 
M”*' and the components of X are redundant coordinates on it. Also, since the left 
S'0(n)-action on Vn,r {X i—)• RX, R E SO{n)) is transitive, it can be realized as a 
homogeneous space SO{n)/SOin — r). 

The cotangent bundle T*Vn,r can be realized as the set of pairs of n x r matrices 
{X,P), P = {pi- ■ -pr) that satisfy the constraints 

X'^X = lr, X^P + P'^X = 0. (3.1) 

The above realization of T*Vn^r is motivated by the natural generalization of the 
Neumann systems on Vn^r given below, however there are other natural realizations 
of T* 14 ,^ (see H]). 

The canonical symplectic structure u) on T*Vn,r is the restriction of the canon¬ 
ical 2-form loq = XlILi del ™ ambient space For our 

purposes it is convenient to work with the redundant variables {X,P), in which the 
corresponding canonical Poisson structure {•,•} on T*Vn^r was explicitly described 
in |23j by using the Dirac construction |151140] . The Hamilton equations have the 
matrix form 

dH ■ dH , , 

^ = P = -— + XA + Pn. (3.2) 

where A and H are the rxr symmetric matrix multipliers uniquely determined from 
the condition that the trajectory {X{t), P{t)) satisfies the constraints (|3.in . 

The Lie groups SO{n) and SO{r) naturally act on T*V{n,r) by left and right 
multiplications respectively, with the equivariant momentum mappings given by 

: T*V{n,r) —>• so(n), = PX^ — XP"'" {= pi A ei + ■■■+ pr A e^), (3.3) 
T : T*V{n, r) -A so(r), T = X^P - P'^X. (3.4) 

Here we identified so(n) = so(n)* and so(r) = so(r)* by the use of the invariant 
metrics on so(n) and so(r) dehned by {r]i,r] 2 ) = 

Note that the actions of SO{n) and SO{r) commute and, in particular, the 
components of the momentum mapping, and 

^ij = {ei,Pj) - {ej,Pi), 1 <i <j <r, (3.5) 

are SO{r) and S'0(n)-invariant functions, respectively. 

The Neumann systems. By analogy with the system ()2.7I) on the sphere 
one dehnes a Neumann on the Stiefel variety Vn,r as a natural mechanical system 
with an S'0(n)-invariant kinetic energy and the quadratic potential function 

T 

P =--tr(A^AA) =--Aej), A = diag(ai,..., a„). (3.6) 

i=l 
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Following [23], we consider a family SO{n) x S'0(r)-invariant metrics defined 
by the kinetic energy 


r,(X, P) = ch) + ^) = i tr(P^P) - 


1 


+ K]iT{{X^Pf), 


K being a parameter (k > — 1). In the class of the metrics ds‘^ we have the normal 
metric induced from a bi-invariant metric on SO{n) {k = 0) and the Euclidean 
metric induced from the Euclidean metric of the ambient space (k = —1/2). 

The corresponding Hamiltonian flows are given by 


X = P-(1 + 2k)XP^X, 

P = AX + {1 + 2k)PX'^P + XA, 

where the Lagrange multiplier matrix does not depend on n: 

A = -X^AX - P^P. 


(3.7) 


(3.8) 


In particular, the Neumann system with the Euclidean metric is given by 

X = P, P = AX + XA. (3.9) 

Note that, due to the S'0(r)-symmetry, both of the potential (13.6p and the 
metrics dsf., the momentum mapping T is the integral of the system (|3.7I) . The 
difference between the two vector fields in (13.7p with different k is proportional to 
X' = -XP'^X, P' = PP^X. Since X^P = -P'^X = the latter describes 
permanent (steady state) rotations of the vectors ei,..., e,. in = span(ei,..., e,.), 
similarly for the vectors pi,... ,Pr. 

The geometry of Riemannian spaces {Vn,r,ds‘^) is studied in [28]. It appears 
that the class of the metrics ds\ is suitable for studying other natural mechanical 
problems on Vn^r such as a variant of the pendulum with a linear potential V = 
( 7 I) Cl) + • • • + (jr, Sr) (ji,... ,jr are constant vectors in M"), as well as the 4th 
degree potential V = ti{X'^A^X) — tT{X'^AXX'^AX) (see [23]f. 


The ”big” (nxn) Lax representation The derivatives of the S'0(n)-momentum 
= PX"^ — XP^ and the symmetric matrix XX"^ along the flows (13.71) are given 
by: 

A], = [‘I>,XX^] . 

As a result, the following theorem holds 

Theorem 3.1 [23] Eor all k, the equations (1321) admit the same Lax representation 

4£(A) = |Af(A),£(A)] (3.10) 

with a spectral parameter X, where C{X) and A7(A) are nxn matrices given by 

£{X) = X^ + XX^+ X‘^A, A7(A) = $ + AA. (3.11) 
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The Lax representation (13.101) . (13.111) is closely related to that for the integrable 
Clebsch-Perelomov rigid body system [l5] and for r = 1 it was given by Moser in 
|4U| . It belongs to the class of the Lax matrix representations related to symmetric 
pairs decompositions of Lie algebras m- 


The ’’small” (2r x 2r) Lax representation. 

Theorem 3.2 [23] Up to the action of the discrete group "Ltf generated by n reflec¬ 
tions 


{e\,..., el,p\, ...,pI) — i-e \,..., -e^, -p\,..., -pi), 


i = l,...,n (3.12) 


the Neumann flows (1331) are equivalent to the following 2r x 2r matrix Lax repre¬ 
sentations with a rational spectral parameter A 




L{X) = 


= [L{X),N^{X)], 

( X^{Xln-A)-^P 
VI. - P^(AI„ - A)-^P 


xT{Xln-A)-^X\ 
-P'^{Xln-A)-^x) ’ 


N^{X) 


({1 + 2k)X^P I. A 

V AI. + A -{1 + 2k)P'^x) ' 


(3.13) 

(3.14) 

(3.15) 


the matrix r x r factor A being already defined in (ITSD . 

The statement is checked straightforwardly by using the constraints (|3.ip and 
the matrix identity 


A(AI„ - A)-^ = (Ain - A)-iA = A(AIn - A)"^ - I„ . (3.16) 


Note that the small Lax pair ()3.13p - ()3.15p is a direct generalization of the 2x2 
Lax pair (|2.8p . (12.9p . and it was first given in unpublished manuscript [33|. Both 
Lax matrices, ()3.14p and ()3.1ip . produce the same hrst integrals and their spectral 
curves are birationally equivalent by the Wienstein-Aronszjn formula (see mm)- 


Integrability. From the Lax matrix T(A) we obtain the set of commuting integrals 
S' = {fk,iiX,P)} defined by 

ti{X{PX^-XP^) + XX^ + X^A)’^ = k = l,...,n. (3.17) 

i 

Apart from these integrals, which are S'0(r)Anvariant, the Neumann flows also 
possess the non-commutative algebra of integrals 21 = T*(]R[so(n)]) of polynomials in 
Tjj. As was shown in [23|, the above systems are integrable in the non-commutative 
sense (i.e., they are superintegrable, [431138]) and the dimension of their invariant 
isotropic tori is less than the dimension of VA,.- 

To describe the geometric structure and the dimension of the invariant tori, 
consider the Poisson bracket II on gl(n), obtained from the canonical Lie-Poisson 
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bracket in the dual to the semi-direct product so(n) 0ad Sym(n), via identification 
of vector spaces gl(n) = so(n) 0 Sym(n) = gl*{n) [23]. Then 

0 : (T*K,r, {, }) ^ (gl(n), n), {X, P) ^ 0(X, P) = PX^ - XP^ + XX^ 

is a Poisson mapping, which is invariant with respect to the above right SO{r)- 
action on {X,P) and the reflection Z 2 : {X,P) 1 —>■ {PX,PP). Thus, 0 can be 
regarded as the factorization 


T*Vn,r e{T*Vn,r) = T*Vn,rI SO{r) 1^2. 


Note that the integrals from the center J of 21, generated by S'0(r)-invariant poly¬ 
nomials tr('k^^), /c = 1,..., [r/2], induce the Casimir functions J^, k = 1,..., [r/2] 
of the Poisson space Q{T*Vn,r)- 

The reduced system on 0(14, r) C gl(?^) can be obtained as a restriction of an 
integrable n-dimensional top with symmetric inertia tensor moving in quadratic 
potential (in particular, the Bogoyavlenski generalization of the Clebsch-Perelomov 
system m ). However, the image Q(T*Vn,r) is the union of singular sympleciic leaves 
in (gl(n),n). Namely, the generic symplectic leaf in (gl(n),n) has the dimension 

— n, while the dimension of a generic leaf in 0(T*14,r) equals 


21 = 2r{n — r) + 


r{r - 1 ) 

2 


T' 

. 2 . 


< — n . 


(3.18) 


Thus, the integrability of the reduced Neumann systems does not follow directly 
from that of the above top. In [23] we proved that if all the eigenvalues of A are 
distinct, then the 50(r)-invariant integrals 5^ induce a complete commutative set on 
a generic 2Tdimensional symplectic leaf Uc = {Jk = Ck} C Q{T*Vn,r)- As a result, 
there are ^0r(r — l)/2 independent functions in 0 21 and its center 0 J has 
I 0 [r/2] independent function^. Since dim T*Vn,r = {I + — l)/2) + {I + [r/2]), 

by the theorem of the non-commutative integrability (see e.g., [3I1I381I13]), the set 
of polynomial 5^0 21 is a complete noncommutative set of functions on T*14,r- This 
leads to 


Theorem 3.3 [23] Let all the eigenvalues of A be distinct. Then the Neumann sys¬ 
tems (ISZD are completely integrable in the non-commutative sense with the commu¬ 
tative set of integrals given by (|3.17p and by the components of the SO{r)-momentum 
mapping ([331). The generic motions of the systems, with the momentum T of the 
maximal rank, are quasi-periodic over the isotropic tori of dimension 


5 




r) 0 


r(r — 1) 
2 



(3.19) 


Note that there is also an alternative path to the reduced system. Namely, 
the symplectic leaves of the Poisson reduced space T*Vn,r/SO{r) are the Marsden- 
Weinstein reduced spaces ^~^{h)/SO{r)h, where SO{r)h = {Q G SO{r) \ AdQ{h) = 

^In fact, it can be proved that the invariants 3 are linear combinations of the polynomials in J. 
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K} is the isotropy subgroup of /i E so(r) (e.g., see [S^)- For a regular h, SO{r)h 
is the maximal torus The factorization of by SO{r)h = and 

by the reflection Z 2 : {X,P) —>■ {—X,—P) coincides with the corresponding leaf 
Uc C Q{T*Vn,r) with Jk = Ck, Ck = , [r/2\. 

In this way we obtain the rednction of the restricted Neumann flows on T“^(/i) 
to the symplectic manifolds Uc- The systems on Uc are integrable in the usual 
commutative sense: the invariant isotropic tori laying on reduce to the l- 

dimensional Lagrangian tori laying on Uc- These observations can be snmmarized 
in the following diagram, where the vertical arrows represent factorizations and the 
horizontal arrows the corresponding inclusions: 


T*V{r, n) 

ISO{r)IZ2 

gl{n) D &{T*V{r,n)) i 


'i>=h 


Jl=Cl, J[r/2]=C[r/2] 




/Th/2] 


Uc D T^ 


We do not consider intermediate cases, when rank of T is not maximal, and 
proceed straight to the case T = 0 . 


The case T = 0 and reduction to the Grassmannian variety. The oriented 
Grassmannian Gn,r is the variety of r-dimensional oriented planes passing through 
the origin in M”. It is a qnotient space of the Stiefel manifold Vn,r with respect to 
the right S'0(r)-action via submersion 7 r(ei 62 • • • e^) = ei A 62 A • • • A e^, and its 
cotangent bundle T*Gn,r is symplectomorphic to the reduced space 'I'“^(0)/5O(r). 
Note that {X,P) belongs to 'I'“^(0) if and only if X^P = 0. Thus, the reductions 
of the Neumann systems have the same kinetic energy given by a normal metric on 
Gn,r for all K. 

The image 0('I'“^(O)) is a single, 2(n — r)r-dimensional symplectic leaf Uq in 
(gl(n),n). Thus, 0 establishes an isomorphism between T*Gn,r /'^2 and WqH If all 
eigenvalues of A are distinct, then the integrals ^ factor to a complete commutative 
set on Uq (see [23]). Therefore, the rednced Neumann flow on T*Gn,r is completely 
integrable in the usual Liouville senseEl 


4 The spectral curve 

The Lax matrix (|3.14p is a particular case of so called rank p perturbations of the 
constant matrix A = diag(ai,..., an), which generate Lax pairs of a great variety 
of integrable systems and have the form 

n .r 

L{\) = Y + Y,Tr^. y,A/*Egi(p), (4.1) 

® Note that the Z 2 -action on Gn,r is trivial for even r: ei A • • • A er — > ( —l)’'ei A ■ ■ ■ ACr- 

can be proved that on the unreduced space 4'~'^(0) the invariant tori of the Neumann flows 
are also r(n — r)-dimensional. That is, over each torus in the factor 4'“^(0)/5O(r) there is SO{r)- 
parametric family of tori in 
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where V is constant and Afi depend on the variables of the corresponding dynamical 
system. In our case p = 2r and 




^f^ = 


-p\ 


{p\ 


Pr 


4) 


Y = 


0 0 
I. 0 


(4.2) 


\-Pr/ 


General cases of rank p perturbations corresponding to Y with a simple spectrum 
and the properties of the corresponding spectral curves were studied in detail in 
[13 ill la [3]. The matrix (j3.14p however has a quite special structure, and its 
spectral curve S has extra strong singnlarities. So, the previous results concerning 
its genns, number of infinite points, etc, should be adopted to our case. For this 
reason below we give our proper self-contained description. 

Namely, multiplying L{X) in ()3.14D by a(A) = (A — oi) • • • ~ cin), we obtain 

the polynomial Lax matrix 

L(A) = a(A)L(A), (4.3) 

defining the spectral curve 


5 C C^{X, w} : F{X, w) = det(L(A) - wl 2 r) = 0, (4.4) 


which is birationally equivalent to the spectral curve of the Lax matrix C: det(/l(A) — 
vln) = 0. By expanding F{X,w) in w we get 

Fix, w) = + u; 2 "- 2 a(A)X 2 (A) + • • • + iX)l 2 r- 2 {X) + a‘^^-\X)l 2 riX), 


where l 2 i{X) is a polynomial of degree n — / in A with the leading coefficient Cf = 
Their explicit expressions are given in [23] FI 
Due to the symplectic block structure of L(A), the coefficients at odd powers of 
w in T(A, w) are zero, so the spectral curve has the involution 


a : (A, w) (A, —w). 


(4.5) 


Note that, although the Lax matrix L(A) is not invariant under the right SO{r)- 
action, the spectral curve and therefore all the integrals F 2 i{X) are invariant. As a 
simple counting shows, the number of nontrivial integrals on T*Vn,r provided by the 
coefficients of I 21 W eqnals (n — 1) -|- (n — 2) (n — r) = r(n — r) + r(r — l)/2, 

which coincides with the dimension of the Stiefel variety. Since the number of 
integrals is bigger than half of the dimension ()3.18p of a generic symplectic leaf 
within {T*Vn^r)/SO{r), some of the integrals are dependent. Like the ’’big” Lax 
matrix C{X) in ()3.10p . the small Lax matrix L{X) does not produce explicitly the 
momenta integrals Tjj, which should be counted apart. 

®One should replace A and A by —A and —A, respectively, to relate the small Lax matrix L{X) 
and integrals X 21 A) used here with the matrix CneumW and integrals X 2 ii\) used in |23 |. 
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4.1 Genus of the regularized spectral curve. 

From now on we consider the compactified version of the curve (I4.4p . identifying it 
with its projective closure in = (^ : ry : C) such that A = w = r]/(. 

First, note that in its finite part, S has singular points Oj over X = ai, i = 
1,..., n, where all the branches of w{X) meet and all the partial derivative of F(X, w) 
vanish up to order 2r — 2. 

To perform a desingularization of S here, we use the fact that the eigenvectors of 
the polynomial matrix L(A) for X = ai are proportional to those of A/i in (14.Ih . As 
was shown in [3], if rank Mi = ki, then L(aj) has 2r — fc, independent eigenvectors. 
In our case, due to (14.2p . all Mi have rank 1 and there are 2r — 1 independent 
eigenvectors. Hence, the regularized curve S' in the projection 

S' ^ P{A} 

has 2r — 1 different points and one ordinary (second order) branch point Pq. over 
A — ai. 

Following the theory of singularities (see e.g., m), this implies that at Oj, the 
curve S has singularity (d-invariant) of order 

5i = (2r - l)(r - 1). 


Singularity at the infinite part. The projective closure of 5 in P^ = ((^ : ry : 
has only one infinite point (0:1:0). Indeed, the structure of F{X,w) in (14.41) 
suggests that w = O (A^^"”^)/^) as A —)• oo. 

Then, in the neighborhood of (0 : 1 : 0), apart from the local coordinates C, it 
is convenient to use the weighted coordinates t, lu such that 

1 tu F 

A = w = — and, therefore, f =-, C = — • (4.6) 

t F ’ ’ ^ tu tu ^ ^ 

Substituting this into (14.41) and multiplying by we obtain the equation 

tu^'’4-1- , - \-{-tY+OriKl,t) = (tU^-t)'’+Or(tU,t) =0, (4.7) 

t!(r — /)! 

where Or(tu,t) denotes higher order terms tu“t^ with 2 q; + /3 > r. Thus the point 
(t = 0, tu = 0) indeed belongs to S. The above equation also shows that S has a 
strong singularity at the infinity. 

To desingularize the curve there, consider the expansion of the polynomial Lax 
matrix dO]) 


L(A) 


(V{X) U{X) \ 

VW(A) -V^(A); 

( A’"-^Vo + A’"-2Vi + • • • X^-^lr + X''-‘^Ui + ■ ■ ■ \ 

VA^I^ + A^-^Wo + A"-2>Vi + • • • -A’^-^V^ - A”-2vf + ■■■) 


(4.8) 
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with the leading r x r matrix coefficients 
Vo = X^P, Vi = X^AP, Ui = -tiAlr + X^AX, Wq = -tr Air - P'^P , 


which are r x r matrix generalizations of the coefficients (I2.13p . Note that the 
components of the coefficient Vq = X^P are first integrals of the system: Vq = 
Next, substituting (14.61) into the eigenvector equation 


L(A)V' = wif:, G 


we get the expansion 


/ V^t + Vit^ + • • • 
VI, + Wot + Wit^ +... 


I,t lA\t^ 4“ * * * \ 

+ ■■■) 


V’(t) = tt)(t)V’(t). 


(4.9) 


Hence, as t —>■ 0, the eigenvectors of L(A) tend to those of Y 



which has 


r different eigenvectors with zero eigenvalue. 

As a result, the regularized curve S' has r different points over A = oo, which 
are all ordinary branch points and will be denoted as ooi,..., oo,. 

The structure of (j4.9p allows to evaluate the expansions tv{t) and in the 
neighborhood of each oog- 


Proposition 4.1 1) Let rank of Vq = X^P ^ 0 be maximal and r is odd. Then 

let 

{z/l, . . . , • • • - W/2] ’ 

be the set of the eigenvalues of the matrix coefficient Vo G so(r) and 


{vi, . . . V[,/ 2 ], Vi, . . . , V[,/ 2 ],Vo} 
be the corresponding eigenvectors normalized by 


Wxs) = l, {P,Vs) = l, (/3,vo) = l, 

where /3 = (/3i,...,/3,) G C*” is any non-zero vector, not orthogonal to any 
eigenvector o/ Vq. Then in the neighborhood o/ooi,... ,oo, the expansions of 
Xo{t) are 

tVs{t) = + Vgt + hgt^^'^ 4-, 

X0]^r/2]+s{t) -^st + bst^/"^ -, s = l,...,[r/2], (4.10) 

tu,(t) = 4- b^t^/"^ + Bot^/^ 4-, 

®We choose the indices of the coefficients in the same way as in (ESP 
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and the corresponding expansions of normalized by 
l3iip[r/2]+i H-K l3r'if[r/2]+r = 1 are 


Ipsit) 


V’[r/2]+s(0 

ifrit) 


( Vs + hgt + 

V Vs + hsi^/^ + 0{t) J ’ 

/vsi^/^ + hsi + 0(i3/2)\ 

V Vs + hsi^/^ + 0{t) J ’ 

/vot^/2 + 0(i3/2)\ 

V Vo + hot^/^ + 0(i) y ■ 


(4.11) 


Here bs,bo,BQ are certain constants, and hs,hs,ho are constant vectors. The 
last expansions tt)r(t), V’r(^) contain only odd powers oft^^'^. 

2) If rank of Vo is maximal and r is even, then item 1) holds without the last 

expansions. 

3) If Vo = 0, then all the expansions contain only odd powers oft^^'^. 


Proof. It is sufficient to substitute the above expansions into (|4.9p and compare the 
few leading terms in both sides. □ 

The involution a : S ^ S acts on the expansions as (t, {t, —to{t)). On 

the other hand, due to the presence of each expansion of Proposition [4T] gives 2 
branches of S' near oos. Then, if rank of Vo = X^P is maximal, for s = 1,..., [r/2], 
the second branch of rDs(t) obtained by replacing y/t ^ —y/t coincides with the 
first branch of —tD[r/ 2 ]+s(^) a-nd vice versa. If r is odd, the two branches of Wr{t) 
transform to each other under the change tu —)• — to. 

As a result, the infinite points of S' admit the natural division into the subsets 

{ooi ) ■ • • o®r/2}’ 

{c>Oi , • • • 

such that a{ooJ) = oo^, j = 1,... , [r/2] and it(ooo) = ooo. 

If Vo = 0, as follows from item (3) of the above proposition, all the infinite points 
ooi,..., oor are invariant with respect to the involution. 

Using the expansions tt)s(t) and the second pair of relations (|4.6p . one calculates 
the set of the Puiseux expansions of the coordinate C in powers of ^ in the neigh¬ 
borhood of the infinite point (0:1:0). Then, applying the theory of singularities 
of algebraic curves (see e.g., mm), we calculate the order 6oo of singularity of S 
at (0 : 1 : ofl: 

6oo = 2nr(nr — 2r — 1) -|- 2r(r -|- 1) if rank (A^P) is maximal, 

6oo = 2nr{nr - 2r - 1) + r"^ + ^ if = 0. 

^These long calculations are due to Maria Alberich (UPC) 


{ooy,...00+2} 


if r is even. 


{ooy,... ooi^/ 2 ]}’ ^ odd. 


(4.12) 
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Now, summing the singularity orders 600 and (5i,..., 5^ and taking into account 
that the degree of the curve S equals N = n{2r — 1) + (n — r), we use the Pliicker 
formula to get the following result. 


Theorem 4.2 The geometric genus of S equals 


{N - 1){N - 2) 
2 


g = 2 r{n — r) — n + 1 


3 1 

- 600 = 2 nr -n - - -r + 1 

i=l 

if rank (X'^P) = 2[r/2], 
ifX^P = 0. 


(4.14) 


4.2 The ’’small” curve C and its genus. 

In view of the involution (|4.5I) , the curve 5 is a 2-fold ramified covering of the curve 
C C C^(u, A), u = , given by the equation 

C : u" + u^-^a{X)X2{\) + ••• + «• a^^-^{X)T2r-2{\) + a^^-\\)T2r{X) = 0 


Making the birational transformation (A, u) (A, /r) with u = a(A) //, we see 
that the equation of C gives the lines {A = a*} and the curve 

C \ fp if ^X2(A) P ■ ■ ■ P fid'' ^(A)X2r.— 2 (A) -|- d’’ ^(A)X2r.(A) = 0. (4.15) 


The latter is singular over A = ai,..., and at its infinite part. Let C be a compete 
regularization of the curve C. Then the regularized curve S' is a 2-fold covering of 

C. 

The covering vr : S' —?■ C is branched when // = 0, that is, when A = a, or A is 
a simple root of the last polynomial X 2 r{X) of degree n — r in (|4.15|1 . As mentioned 
above, the projection 5' —)• P = {A} has only one ordinary branch point Pa^ over 
A = Oj, hence Pa^ is also the only branch point of the covering vr over A = a*. (This 
means that the projection C' ^ = {A} is not ramified over A = a*.) The roots of 

X 2 r{X) give other n — r ordinary branch points of vr. Therefore, tt has always 2n — r 
finite ordinary branch points. 

Next, as follows from the expansions to{t) of Proposition 14.11 if rank('k) = 2[r/2] 
and r is odd, the covering tt is also ramified at only one of the infinite points 
ooi,..., oor, where the two branches of tOs{t) pass to each other under the involution 
a, and it is not ramified over infinity when r is even. 

If 4' = 0, then, following item 3) of Proposition 14.11 the two branches of tOs{t) 
pass to each other under a for any s = l,...,r, which gives r infinite ordinary 
branch points of tt. 

As a result, the covering vr : S' —)• C has in total 2 n — 2[r/2] ordinary branch 
points if rank(4') = 2[r/2], and 2n ordinary branch points if 'L = 0. Then, according 
to the Riemann-Hurwitz formula, 


2rj — 2fr-/2l 

g = gen(5') = 2(gen (C) - 1) +-+ 1, if rank(4') = 2[r/2], 


5 = 2(gen (C')-l) + —+ 1, 


(4.16) 


if 4- = 0. 
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In view of (14.141) . this gives, 


Proposition 4.3 If rank{^) = 2[r/2], 
90 = gen{C') = 

and, i/ 'h = 0 , 


= 2 [r/ 2 ], 


J n(r — 1 ) — 3(r^ — l)/4 r is odd, 

1 n{r — 1 ) — 3r^jl: + 1 r is even 

(4.17) 

go = r{n — r) + n — 1 . 

(4.18) 


In the next section the genera of S',C' will be used to calculate the dimension of 
the Abelian subvarieties of the Jacobian of S', which are related to complex invariant 
tori of the system. 


The case r = 2. In this simplest case it is possible to calculate the genera gQ,g 
without tedious calculations based on the above Puiseux expansions. Namely, the 
curve C in (I4.15P takes the simple form 

C + /iX2(A) + a(A)X4(A) = 0. 

which, under the birational change g, = —X 2 (A )/2 + y/2, gives the Weierstrass 
hyperelliptic form 

=Xl{\)-Aa{X)X^{X). (4.19) 

To find its genus, note that here X 2 (A),Zi(A) are polynomials of degrees n — 1, n — 2 
respectively, hence one might expect the right hand side of (|4.19p to be a polynomial 
of degree 2n — 2. However, in the general case (when 'k 7 ^ 0), the degree is 2n — 3 = 
2 (n — 2 ) + 1 (in view of the leading coefficients of Z 4 (A),X 2 (A), the coefficient at the 
leading power vanishes). Then the regularization C' of C has genus go = n —2 

and one infinite point. 

Next, as was counted above, for r = 2 and 'k 7 ^ 0 , the covering vr : S' ^ C 
has in total 2n — 2 ordinary ramification points. Hence, by the Riemann-Hurwitz 
formula (j4.16p . the genus of the ’’big” curve S' is 

g = 2(n - 2 - 1) + (2n - 2)/2 + 1 = 3n - 6 , 

and S' has 2 infinite points over A = 00 , which pass to each other under the involution 

a. 

In the case 'k = 0 the right hand side of (j4.19p becomes a polynomial of degree 
2 n — 4, the curve C has now 2 infinite points, its genus drops to go — 3, whereas the 
covering vr has 2n ordinary ramifications. All this yields g = 3n — 7. 

Remark 4 In the general case r > 2, the curve C is not hyperelliptic, and to 
calculate its genus directly (i.e. without relating it to the upper lying curve 5), we 
need to know the number of finite points of ramification of C' —>• P = {A}. This 
could be done by calculating the degree of the discriminant of the polynomial (j4.I5p 
in g, however, due to quite special values of the leading coefficients of X 2 i{X) in (14.411 . 
the correct degree is much less than the expected one, and is quite hard to be found 
for the general n,r. This effect is due to the strong singularities of 5,C at their 
infinite points. 


21 














5 Complex invariant manifolds 


From the curve S' we pass to its Jacobian variety, Jac(5'), defined as the additive 
group of degree zero divisors on S' modulo divisors of meromorphic functions. Equiv¬ 
alently, in some cases we will consider effective divisors D as elements of Jac(5') by 
choosing a basepoint Pq ^ S' and associating to D the degree zero divisor D — NPq, 
where N is the degree of D. 


5.1 The Prym variety 

The involution a : S' ^ S', a{X, w) —)• (A, —w) extends to the Jacobian of S'. Hence 
Jac(5') contains two Abelian subvarieties: the Jacobian of the ’’small” underlying 
curve C of dimension go= genus C and Prym(5', a) of dimension k = gen(5') — 
gen(C'). The vectors of Jac(C') C Jac(5') are invariant with respect to a, whereas 
the vectors of Prym(5', a) C Jac(5') are anti-invariant. 

Algebraically, Prym(5',cr) is defined as the set of degree zero divisors P on S' 
such that 'D + a'D is the divisor of a meromorphic function on S'. For a given degree 
zero divisor P, a translate of Prym(5', a) in Jac(5') is the set of degree zero divisors 
V satisfying 

V + aV = n, 

where = denotes the equivalence modulo divisors of meromorphic functions on S'. 

The variety Jac(5') is isogeneous (but not isomorphic) to Jac(C') x Prym(5', a). 
A detailed algebraic description of Prym varieties in the case of a general double 
covering S' —)• C ramified at N > 2 points is given in [19] (Chapter V). Following 
this description and our previous calculations, Prym(5', a) is an Abelian variety 
with polarization (!,...,!, 2 ..., 2 ) if rank T is maximal and (1,... , 1, 2 ..., 2) if 

go n-[r/2]-l go n-l 

I' = 0. 

From the expressions (I4.14|] . (j4.17p . and (I4.18|) we immediately obtain 

Proposition 5.1 1) For generic constants of motion, the dimension k of the 

Abelian subvariety VT:ym.{S', a) equals 



-r) + 


r{r - 1) 
2 



which coincides with the number of independent and commuting integrals ^ and 
with one half of the dimension of a generic symplectic leaflAc in &(T*Vn,r)- 

2) In the special case 'k = 2X^P = 0 the dimension o/Prym(5', u) is r{n — r), 
which coincides with the dimension of Gn,r- 


Note that in the simplest case r = 2 (and only in this case), the dimension of 
Prym(5', cr) equals 2(n —2), i.e., the dimension of the Grassmanian Gn, 2 , regardless 
to whether 'k = 0 or not. 
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5.2 The generalized Jacobian and afRne Prym variety. 

Assuming that the rank of the momentum 'I' is maximal, we also introduce the 
generalized Jacobian of the singularized curve S" obtained from S' by gluing pairwise 
the involutive infinite points ooj, oof, s = 1,..., [r/2] described in (14.121) . Namely 

S = 5o U {ooi, . . . , 00 |-^/ 2 ]}) = 5 \ {oO]^ , OO ^, ■ ■ ■ , °®jr/ 2 ]} 

SO that odi,... ,oo[r/ 2 ] are ordinary double points of S". Then a meromorphic func¬ 
tion / on S' will be also meromorphic on S" if f{ooj) = f{oof). 

The generalized Jacobian of S", denoted as Jac(5', oo) is then defined as the set 
of degree zero divisors V on the affine part 5o considered modulo oo-equivalence: 
V =oo Q if there exists a meromorphic function f(p) on the regular curve S' such 
that 


l)(/) = -P-Q, 2)/(ooj =/(oo+) 7 ^ 0,oo, s = 1,..., [r/2] 

(note that the values f{oof), /(oo^) for i j can be different)!! 

Thus, there is an exact sequence of algebraic groups 

0 J^(5',oo) Jac(5') -^ 0 (5-1) 

where, for p = {pi,... ,p\r/ 2 ]) € the image v{p) is the divisor of any mero¬ 

morphic function f{p) on S' satisfying 

f{oot)/f{oo~) = Ps, s = l,...,[r/2], (5.2) 

and for any degree zero divisor V, <()(P) gives a point in Jac(5'). 

The analytical description of Jac(5',oo) is following: let uJi,... ,u}g he g inde¬ 
pendent holomorphic differentials on S' and f2j, j = 1,..., [r/2] be meromorphic 
differential of the 3rd kind having a pair of simple poles at oo^,ooJ respectively. 
One can always normalize the meromorphic differentials in such a way that they 
will be anti-invariant under the involution: = —h2j. 

Next, let A C C^{zi,..., Zg) be the lattice generated by vectors of periods of 
(wi,... ,ojg)'^ with respect to a basis in Hi{S','Z). Respectively, let 
A' C ... ,Zg,Zi,..., Zy^i 2 \) be the lattice generated by vectors of periods 

of 

{iOl,...,LOg,Ql,..., 17[^/2])'^ 

with respect to a basis in Hi{Sq,7j) formed by 2 g canonical cycles on S' and the 
homology zero cycles embracing oo/",..., oo[];y 2 ]- The lattice A' has rank 2g + [r/2], 

and the generalized Jacobian Jac(5',oo) is the factor which is a non¬ 

compact algebraic group. Topologically, it is the product of the customary Jacobian 
Jac(5') = C^/A and of [r/2] cylinders C* = C \ {0}, however algebraically its 

®The above definition should be distinguished from the other kind of generalized Jacobians, 
when several points of a smooth curve are glued to the same point, as described in [1011521 [26] . 
Both kinds include the simplest case when the singularized curve has only one double point. 
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structure is more complicated in the sense that its period lattice A' is not a direct 
product of the period lattices of Jac(5') and of The above map (j) is just 

the projection 

= {zi,...,Zg). (5.3) 

The relation between the above analytical and algebraic descriptions of Jac(5', oo) 
is given by the generalized Abel map with a basepoint Pq £ 5*0 

fP 

A{p)= PgSo. 

JPo 

A general point in Jac(5', oo) is the Abel image of an effective divisor of degree 
> g + \r /2] on So ■ The inversion of the generalized Abel map in terms of generalized 
theta-functions is described in [T3l [20] . 

A generalization of the Abel theorem says that two effective divisors 
P = Pi + P/v, Q = Qi + ■ ■ ■ + Qn on So are oo-equivalent (P =oo Q) if and 
only if 

A(Pi) -l- • • • -l- A{Pn) = A{Qi) -!-•••+ A{Qn) 

modulo the period lattice A'. 

Note that due to this definition, two divisors P^P 2 on So can correspond to the 
same point on Jac(5), but to different points on Jac(5',oo). 


The affine Prym variety. The involution (14.511 extends to Jac(5',oo), which 
then contains two algebraic subgroups: one is invariant with respect to a and the 
other is anti-invariant. The first is just the usual Jacobian Jac(C^) of dimension ^o- 
The second is the ajfine Prym suhvariet^ Prym(5',cT) C Jac(5',oo), which can be 
defined as the set of degree zero divisors P on So such that 

P -|- (tP = {a divisor of a meromorphic function f{p) 

with /(oo+) = f{ooj) / 0,oo, s = l,..., [r/2] } 

= the origin in Jac(5^,oo) . 

In view of the definition of the map v : i—)• Jac(5',oo), the action of a 

on the set v ((C*)['’/^l) is given by 

a{v{pi,. . ■,p[r/2])) = vCi^/pl,- ■ ■ , ^/p[r/2]) ■ 


Then we obtain 

Lemma 5.2 The set V ((C*)l'’/^1) belongs to Prym(5' , fj). The involution a on it 
has fixed points u(±l,... ,±1), which are also half-periods in Jac(5',oo) and 
in Prym(5', a). 

^We do not use term ’’generalized Prym variety”, quite natural in this situation, since it is 
related to another construction considered in [9]. 
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Proof. Let V be the divisor of any meromorphic function f{p) on S' satisfying 
dS^I), and let f^{p) = f{a{p)) so that /^(oo+)//^(oo7) = l/ps. Then (/ • /^) = 
V + aV and f{oof)fa-{oof) = f{oof)fa-{ooj), hence, by the definition, V belongs 
to Pryni(5', a). 

Next, if V is the divisor of a meromorphic function f{p) satisfying f(oof) = 
±f{ooj), then 2(/^) = V and /^(oo^) = f'^{oof). That is, 2P corresponds to the 
origin in Jac(S",oo), and P is a half-period of the generalized Jacobian and of the 
affine Prym subvariety. □ 

By the above lemma, the exact sequence (I5.ip implies the following sequence of 
groups 

0 {C*)^r/ 2 ] P^{S',a) Prym(5,cj) -^ 0. 

Since (f{v ((C*)!'"/^!)) = 0, the sequence is also exact. 

It follows that, like Jac(5',oo), the variety Prym(5',cj) is non-compact and has 
dimension g + [r/2\ — go = Z -|- [r/2], which, in view of Proposition 15.11 equals 6 , the 
dimension of generic invariant tori in T*Vn,r- 

The generalized Jacobian Jac(5', oo) is isogeneous (but not isomorphic) to the 
product Jac(C) x Prym(5', cr). 

5.3 The eigenvector map 

To proceed further, we need to recall some basic notions relating Lax matrices and 
Abelian varieties, and we follow the description of dziiMiEiiiniEi]. 

Let L(A) he a d X d polynomial Lax matrix and S be its spectral curve (with 
its infinite points), which is assumed to be regular or appropriately regularized, of 
geometric genus g. L(A) defines the eigenvector bundle e : S —>• given by the 

eigenvectors tpiP) = (V'^(P), • • • ,V’^(T’))^, P G S of L(A). 

Impose a normalization (a, '4>{P)) = 1, where a = (ai,..., a 2 rY' ^ is an 

arbitrary constant vector, and consider the minimal effective divisor on S such 
that 

^V’^(P)^ = zeros of V'^(P’) - poles of V’^(P) > —Va, I = 1 ... ,d, 

that is, each of the normalized components fj^P) can have poles at most at the 
points of 'Da- These points can be found as common zeros of certain polynomials of 
A, w (see, e.g., my Also 

deg(Pa) = g + d- 1 = N. (5.4) 

(One can always chose such a normalization a that d — 1 points of Da will be fixed 
in the infinite part of S.) 

It is known that two divisors Da, Da' corresponding to different normalizations 
a, a' are linearly equivalent. Therefore, for a basepoint Pq G S, the degree zero 
divisors Da — NPo,Da' — NPq give the same point in Jac(5'), more precisely, in the 
open subset Jac(5) \ 0, with 0 € Jac(S') being a translate of the theta-divisor, the 
Abel image of all the special divisors on S. 
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Let now Xg be the isospectral manifold: the set of all the above matrices L(A) 
having the same spectral curve S. Thus we get the eigenvector map £ : T 5 —)• 
Jac(5). 

If the equivalence class {P} is the image of a matrix L(A) G X 5 , the latter can 
be reconstructed up to a conjugation by an element of the group PGL(d, C) (not 
depending of A). The main steps of this are as follows: Let L{'D) be the vector 
space of meromorphic functions f{P) on S with (/) > P (it includes / = 1). 
According to the Riemann-Roch theorem, for a non-special degree N divisor P, 
dim (L(P)) = N — g + 1, hence, in the considered case, dim (L(P)) is precisely d. 

There exists a basis {/i(P),..., /d(P)} in L{V) such that f(P) = (/i,..., fdf^i 
P = (A,^) G 5 is an eigenvectoi0 of L(A). To reconstruct L(A), choose A G C with 
distinct eigenvalues /ri,..., and consider the d x d matrix 

F(A) = (f(Pi)---f(P,)), P, = {X,hj)- 

Then the matrix 

A’(A) = F(A) diag(/ii, F-i(A) (5.5) 

is independent of the order of /ri ,..., fJ-d, is, in fact, polynomial in A, and has the 
prescribed spectral curve S. It reconstructs the matrix L(A) up to a conjugation by 
a constant matrix (see e.g., [IZ112]). 

Clearly, two matrices L(A),L(A) G P 5 that are conjugated by a constant element 
of PGL((i, C) have equivalent divisors P,P. The inverse is also true: P = P means 
the existence of a meromorphic function g{P) on S with (g) = P — P. Let the 
components of f(P),f(P) form bases of L(T>), L{'D) respectively. Then the com¬ 
ponents of g{P)i{P) form a basis of L(P), and there is a non-degenerate matrix 
R G GL{d,C) independent of P such that f(P) = g{P)Ri{P). 

It follows that the induced map 

Ai : {IL(A) G Xs up to conjugation by a matrix of PGL((i, C)} 1 —)• Jac(S') \ 0 
is injective. 

In the general important case when the leading coefficient of L(A) is a constant 
diagonalizable matrix J and there are no constraints on the other coefficients, the 
subgroup is the stabilizer of J, i.e., the product (C*)'^”^. Then there exists the 
factor variety Xs/ (G*)^~^ and, as was shown in many publications, if the curve is 
smooth, the induced eigenvector map 

Xs/iC*f~^ ^ Jac(S)\0 

is an isomorphism ([laiiniisiiis]). The inverse map is described explicitly by means 
of theta-functions associated to the curve S. 

If the leading coefficient J in not diagonalizable and/or L(A) belongs to a sub¬ 
group of gl{d,C), the above result should be refined (one must consider the equiva¬ 
lence by conjugations by elements of a subgroup of PGL((i, C)). Some cases of that 
have been considered in [271 ESI [I6]. 

is always supposed that f is normalized, i.e., its components do no have a common zero 
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In addition, some general properties of the eigenvector map in our case, when 
L(A) belongs to the loop subalgebra sp(2r, C), were already studied in [2] (Section 
5, case ii). 

5.4 The complex manifolds X, 

Consider the Neumann system on the complexified cotangent bundle T^Vn,r defined 
by the constrains (13.111 in the complex domain. Let us fix generic values of all the 
I + [r/2] independent commuting integrals 

0*{Pi, • • • ,Pi}, tr(4'2)^..., tr(4'2h/2]) 

described in Section [3l This also fixes the spectral curve S' of the ’’small” 2r x 2r 
Lax matrix (I3.14p . (14.3p . which belongs to the loop subalgebra sp(2r, C). 

Let X C T^Vn^r be the corresponding isospectral manifold, which then has the 
(complex) dimension 


6 + r(r — l)/2 — [r/2] = I + dim SO{r) 


with 6 given by (13.191) . 

The manifold X and the curve S' are invariant with respect to the right action 
of the complex group 50(r, C) on {X,P), as well as to the action of the group Z 2 
of reflections (I3.12p . Moreover, the product SO{r,C) x Z 2 is the maximal group 
leaving X invariant. The reduced manifold = XjSO{r,C )is the complex 
extension of the corresponding /-dimensional invariant isotropic torus of the reduced 
Neumann system on T*V{n^r)jSO{r)= P>{T*V{n^r))■ 

Next, the complex isospectral manifold X itself is foliated by complexihed 5~ 
dimensional invariant isotropic tori 


Xft =xnT-^(/i) = {(X,P) gX I T = /i}. 


the joint level varieties of the extra non-commuting integrals of the set 

As follows from the above, the factor of each X^ by the direct product C* x • — x C* 

[r/2] copies 

which is the complex stabilizer of the momentum h G so(r, C) in SO{r, C) and by 
the group of reflections Z 2 coincides with Xi-ed- This can also be seen as follows: the 
action 


yields 


X^Xg, P^Pg, g£ SO{r, C) 



L(A) = 


°T ) L(A) 


f A^-iVo + --- A’*-il^ + --- A 

VA^I, + A’^-^Wo + • • • + • • •/ ^ 

+ • • • X^-X + ■■■ 

g) V^-I, + A-VWo(7 + --- -A”-VVo^5 + 


(5.6) 


the simplest non-trivial case r = 2 the manifold Th coincides with T. 
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which changes neither the structure of L(A), nor the spectral curve, but may change 
the block Vo = of the so(r)-momenta integrals. If, instead, g belongs to the 
stabilizer of Vq (that is, the latter is preserved as well. Thus the rela¬ 

tions between the manifolds can be described by the following diagram showing the 
factorizations and an inclusion 



A point of 2red = TjSO{r,<C) 11,2 defines an equivalence class {V} of effective 
divisors on S' via the sequence {X,P) i—)■ L(A) V'(T’) i—s- P (where, as above, 
{a,'ll)) = 1). Thus we have the injective eigenvector map 

M : Trod —Jac(5'). 

Note that, in view of Proposition 15.11 the dimensions of Xred and Prym(5', cr) 
coincides. Moreover, we shall prove that M maps the invariant manifold X^ed to the 
translated Prym variety. 


Proposition 5.3 For any equivalence class {(A,P)} in Fred; t/ie divisor T> = 
Ad({(A,P)}) satisfies the relation 


V + aV = B, (5.7) 

where B is the set of all branch points of the projection 5' —^ P = {A} Uoo, including 
the infinite points ooi,, oo^. 

Then, since B is independent of (X, P) G Fred, the degree zero divisor V — NPq 
with a fixed basepoint Pq G S' belongs to a translated Prym variety Prym(5',cr). 
That is, M maps the reduced manifold X^ed into the Prym subvariety. 

Since M is injective and both varieties have the same dimension I, we proved 

Proposition 5.4 The reduced isospectral manifold Fred = X/50(r, Cj/Z^ of the 
Neumann system is isomorphic to an open subset PrymQ(5^,(T) of PTym{S', a). 

Remark 5 In the special case T = 0 the corresponding isospectral manifold X 
has dimension r(n — r) -|- r(r — l)/2, and the factor Fred = FjSO{r,C)fllf is the 
complexification of the r(n — r)-dimensional invariant tori of the reduced Neumann 
system on T*Gn,r- As item (2) of Proposition l5.1l savs. in this case the corresponding 
variety Prym(5^, a) has the same dimension. Then Proposition 15.31 implies that, like 
in the generalcase, the map Ai realized a bijection between Fred an open subset 
of Prym(5', a). 
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Proof of Provosition \5.3[ Here we use the technique already applied in [TTl HU ITT] . 
Namely, let 

f;{P) = ix{P),aP)f, X,eGC^ P={X,w)gS' 

be an eigenvector of L(A) with the eigenvalue w, and ^(P) = ^^{aP) = be 

an eigenvector with the opposite eigenvalue: L(A)'!/j(P) = —'wip{P). Introduce also 
the vector 'if*{P) = x)^, which, due to the structure of L(A) in ()4.8I) . is an 

eigenvector of L^(A) with the eigenvalue w: 

l7{\)f:*{P) =wf:*{P). (5.8) 

Indeed, 

(t -T) (fl=(XI-*;?)=(X) ■ 

Now consider the function 

F{p) = {riP),HP)) = V{p)r{p) = xi- xl ( 5 . 9 ) 

The rest of the proof is based on item 1) of the following lemma. 

Lemma 5.5 1) The function F{P) has simple zeros only at the branch points of 

w as the function of X and has the poles only at V and aV. 

2) F{P) is antisymmetric with respect to the involution a : S' s- S'. 

Since F{P) is meromorphic on S', the divisors of its zeros and poles are equiva¬ 
lent, which leads to the relation (|5.7p and Proposition 15.31 □ 

Proof of Lemma [531 1) Take two distinct eigenvalues w,w' of L(A) and the corre¬ 
sponding eigenvectors if{X,w),if*{X,w'). Then, on the one hand, 

{fj*{X,w')f'L{X)ip{X,w) = w {'ip*{X,w'))^ 'if{X,w). 

On the other hand, in view of (I5.8p . the same product equals 

['ilj*{X,w')Y'L{X)'if{X,w) = {j7'{X)'if*{X,w')f''4){X,w) = w' {'ip*{X,w'))'^'ip{X,w). 

Since w' ^ w, the difference of the above two equations yields {'ijj*{X, w'))'^ ipiX, w) = 
0 for any {X,w) / {X,w') and, by continuity, {ip*{X,w))'^ ^f{X,w) = 0 when {X,w) € 
B. Hence, F{P) vanishes at all the branch points B. 

Next, since the components x(P),.^(P) of tpiP) have poles only at V and, by 
their definition, the components of 'ip*{P) have poles only at aV, the function F{P) 
has poles only at P -|- aT>. Note that the latter points are all finite, that is, distinct 
from ooi,...,oOr G B, if we impose an appropriate normalization, for example 
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Finally, note that, apart from simple zeros at the function F{P) cannot have 
other zeros on 5', because, due to the Riemann-Hurwitz formula, the genus of the 
curve S", as a 2r-fold covering of P, equals g = 2r(0 — 1) + (deg B)/2 + 1, then we 
have deg(^) = 2{g + 2r — 1), and, in view of (|5.4p . 

deg(P + aV) = 2deg(P) = 2{g + 2r — 1) = deg(5). 

2) Since ipiaP) = ipiP) = (x, 0"^ and 'ili*{aP) satisfies L?"{aP) = —w'ip*{aP), 
we have 'ip*{aP) = As a result, 

F{aP) = {r{^P)fi^{<rP) = -xC + xC = -F{P), 

which completes the proof of the lemma. □ 

Since for an appropriate normalization of the eigenvector 'ijj{P) the pole divisor 
V is finite, the eigenvector map M. : Fred —> Jac(5') can be extended to the map 

M : Ih/F^ ^ J^(5',oo) 

Proposition 5.6 The map M is injective. 

Proof. The idea is borrowed from [26]. Namely, let il^{P), P ^ S' he the eigenvector 
bundle of a Lax matrix L(A) G Xh/TTf with the normalization + ■ ■ ■ + = 1, 

and P be the corresponding effective divisor on S' defining a point P — NPq in 
Jac(5', oo). We will show that, given the oo-equivalence class of P, the matrix L(A) 
can be reconstructed uniquely. 

Let D be another effective divisor giving the same point in Jac(5^,oo). Then, 
by the generalized Abel theorem, there exists a meromorphic function G{P) on S' 
such that 

(G)=P-P, G(oo7) = G(oo+) 7 ^ 0,oo, s = 1,..., [r/2]. 

The components of ^^{P) form a basis of L(pj^. Respectively, let ^(P) be a vector 
bundle whose components form a basis of L(P). 

Then the components of G{P)'ip{P) form a basis of L(P), and there is a non¬ 
degenerate matrix R G GL{2r,C) independent of P such that tpiP) = G{P)R'ijj{P). 
The bundles 'ip{P),ip{P) define the same matrix L(A) up to conjugation by R. By 

(IMI), the conjugation preserves L(A) G iff P = , where g G SO{r,C) 

is a stabilizer of Vq G so(r, C). 

On the other hand, the above relation between ip, ip implies 

ip{ook) = G{ook)Rip{ook), k = l,...,r. 

Due to the chosen normalization, near ook the components of ip{ook),ip{ook) are 
finite and given by the expansions (|4.11l) involving the eigenvectors v^, W) vo of Vq. 
This implies 

Vs\\G{oos)gvs, Vs\\G{oos)gvs, s = 1 ,..., [r/2], vo|| G(oor)fl 'vq 
^^Note that these components, in general, take different values at ooj and ooj". 
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hence all vq must be also eigenvectors of g. Since g E SO{r,C), we conclude 

that g and R are the unit matrices, therefore the divisors correspond to the 
same matrix L(A) E X/j/Z^ . □ 

Proposition 5.7 The image AiiXhl'E^) belongs to a translation of Piyni{S', a) in 
Jac(5'). 

Proof. The zeros and poles of the meromorphic function F{P) given by (15.91) imply 
the equivalence (15.711 on Jac(5'), namely V + aV = B, which cannot be extended to 
the equivalence 

2? + (tP =oo B 

because F(oo^) = 0 and the divisor B contains all the infinite points: 

B = Bo + ooi + • • • + oor- This can be repaired by introducing the function G{P) = 
F‘^{P)X, where A is the coordinate on the curve S' C P^(A, tc) and 


(A) — O — 2(ooi + • • • + oOr), 


O being the preimage of A = 0 on S'. Then 

(G) = 2B- 2{V + aV) + 0- 2(ooi + • • • + oo,.) = 2Bo + O - 2(P + aV). 

Moreover, in view of Lemma 15.51 and the symmetry A(crP) = A(P), 

G(oo7) = G(oo+) / 0,oo, s = l,...,[r/2 ]. 

Hence the divisor (G) implies 

2{V + aV) =oo 2Bo + 0 = 2B. (5.10) 

It follows that F — NPq belongs to a certain translate T of Prym(5', a) in Jac(5', oo) 
defined up to translations by half-periods V of Jac(5', oo). 

On the other hand, comparing (I5.10p with (15.71) . one sees any such half-period 
must satisfy </>(H) = 0. Therefore, V includes only the half-periods u(±l,... ,±1), 
which, according to Lemma r5.2l are half-periods in Prym(5', a). Hence, the translate 
of T by P is again T, which proves the proposition. □ 

Since Zh and Prym(5',(T) have the same dimension 5 = 1 + [r/2] and Ai : 
Xh/T+f —)• Jac(5^oo) is injective, we arrive at the main theorem of the section. 

Theorem 5.8 If the rank of the momentum T is maximal, then, after factorization 
by the reflection group Z+f, a generic 5-dimensional complex invariant manifold X^ 
of the Neumann system on Tf.Vn,r is PrymQ(5', a), an open subset of a translate of 
the affine Prym variety Prym(5',cj) C Jac(5',oo). 
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This theorem together with Proposition 15.41 can be summarized with the follow¬ 
ing commutative diagram 


^red 


M 


> Prymo(5',cj) 


M 


> Prymo(5',cT) 


(5.11) 


Remark 6 As was shown in Section 5, in the special case T = 0 all r infinite points 
of the regularized curve S' are invariant with respect to the involution cr. Such a 
situation is described in [9]: one can introduce the generalized Jacobian Jac(5',oo) 
of the singularized curve, obtained from S' by gluing ooi,..., oo^ to one point. In 
contrast to the general case, now the part of Jac(5', oo), which is anti-invariant with 
respect to cr, is just the compact subvariety Prym(5', cr). 


A linearization of the Neumann flows on Prym(5^ cr) and Prym(5', u) will be 
presented in a separate publication. In the sequel we proceed directly to their exact 
discretizations. 


6 Integrable discretization of the Neumann system on 

6.1 The map and its basic properties 

Using the notation of Section [3l consider the following matrix generalization of the 
map (j2.2ip . (j2.16p (discrete Neumann system) with a real parameter A* 

P= Ai/2(A,)A-Ar(A,), 

P = -A^I‘^{K)XpXV{\:), 

where, as above, A(A) = AI^ — A and P is a symmetric r x r matrix, which is found 
from the condition X^P + P"^X = 0 ; 

r(A*) = ^ + X^A^/\XPX^ . (6.2) 

The hrst matrix equation in (| 6 .ip gives the definition of the discrete momentum P, 
whereas the second equation provides the discrete dynamics in an implicit form. In 
view of (16.21) . the structure of (16. ip ensures preservation of the SO{r) momentum: 
JP^P = X'^P. 

Equations (16.ip . (| 6 . 2 p dehne a a multi-valued map (a correspondence) 

: T*Vn,r ^ T*Vn,r, P) = (X, P), 

which can be regarded as a discrete Neumann system on T*I4, J^. 

In fact, it can be considered without the restriction P + X — 0, i.e., can be extended to 
the map on Vn,r{X) x which always preserves the matrix X^ P. 
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(6.3) 


To evaluate the map, we rewrite (16.11) in the form 

1 = A-i/ 2 (A,)(p + xr(A*)), 

P = {P + xr(A*)) r(A*). 

Then, applying the condition X = 1^, we obtain the following quadratic 
matrix equation for T, which defines it implicitly as a function of X, P and which 
generalizes the scalar equation (I2.22h . 

TUT + TV + V^T - W = 0, (6.4) 

where 

U = X^A-^{K)X, Y = X^A-\X^)P, W = lr- P^A-\X^)P. (6.5) 

This is a system of r (r+l)/2 scalar quadratic equations for the r (r+l)/2 components 
of T. However, due to the structure of the matrix coefficients, the number of its 
solutions is less than that predicted by the Besout theorem. The matrix equation 
(j6.4h is known in the literature in connection with stationary solutions of the matrix 
Riccati differential equation, optimum automatic control theory, and its complete 
solution was presented in m- 

In order to describe it, we first observe that the coefficients (j6.5p coincide with 
the r X r blocks of the Lax matrix T(A*) in (j3.14p . Since the latter is symplectic, its 
eigenvalues wi,..., W 2 r are divided into r pairs {wi, —Wi). Let V'l, • • •, V'r G be 
the eigenvectors of T(A*) with distinct eigenvalues wi, .. . ,Wr such that Wi / —Wj- 
The corresponding matrix T = {ipi ■ ■ ■ tpr) will be called a non-special eigenmatrix. 
The following proposition is a direct consequence of the results of Potter m- 

Proposition 6.1 Any symmetric solution of the matrix quadratie equation B) 
has the form 

T ^ 

where H, T are upper and lower r x r halves of a non-speeial eigenmatrix T = 
of L{X^). 



Remark 7 The original paper m solves the equation 

TilT + TTJ + Ti^T - 2T = 0, (6.6) 

where it, QJ, 211 are arbitrary nxn complex matrices. As above, consider the matrix 

^-\fm • 

According to if K has a diagonal Jordan canonical form with distinct eigen¬ 
values wi,... ,Wr such that Wi / —Wj, then the hermitian solutions (f'^ = T) of 
(j6.6p are exactly those described in Proposition 16.11 P = where H,T are 

upper and lower r x r halves of the corresponding eigenmatrix of K. Proposition 
16. II dealing with symmetric matrices can be easily proved following the lines of m- 
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Obviously, the product T is invariant under any gouge transformation 'h i—)• 
^ G, G E GL{r,C). As was also shown in [17|, if 'h is special (that is, contains 
eigenvectors with Wi = —Wj) or its columns are linear combinations of more 

than r eigenvectors ipi, then the product also satisfies (|6.4p . but it is not a 

symmetric matrix. 

Since for generic hnite A* 7 ^ ai there are 2 ^ possible partitions 


W = {rCi, . . . ,Wr I —Wi, . . . , —Wr} 


we conclude that the equation (|6.4p has precisely 2^ complex solutions. They admit a 
natural decomposition {T} = {r_}n{r+} such that for any solution T* E {r_} cor¬ 
responding to a partition w_ = {rci, . .. ,Wr \ —wi, . .., —Wr} there is a unique T** E 
{r+} corresponding to the ’’opposite” partition w+ = {—rci,..., —Wr \ wi, ..., Wr}- 

Since the spectral curve S' has ordinary branch points (A = a*, rc = 0 ), in the 
special case A* = a* there are only r — 1 distinct non-zero eigenvalues wi,, Wr-i 
such that Wi 7 ^ —Wji and the number of symmetric solutions of (16.41) drops to 2”“^. 

Once a solution for T is chosen, X,P are found uniquely from () 6 . 1 h or (j 6 . 8 p . 
Since, for generic A*, different solutions T lead to different images X, P, we conclude 
that the complex map 53^ is 2 ”-valued. 

Continuous limit. Like in the case r = 1, the continuous limit of the map is 
obtained by letting A* —)• 00 . Namely, we set A* = 1/r^ and consider the expansions 


A ^(A*) = (In + -h 0(r^)) , 

X = X + tX + OiT^), P = P + tP + 0{t^). 


This gives 



Indeed, substituting the above into the matrix quadratic equation (|6.4I) and using 
the expansions 


U = T^{Ir + T^X^AX + 0 (t^}), 

V = T^(X^P + t^X^AP + 0(t^)), 

W = Ir- T^iP'^P + t'^P'^AP + 0{t^)), 


we hnd that the coefficients at r*^,r, vanish, which justifies (|6.7p . 
Next, substituting the latter into the equations ()6.3I) . and using 



we conclude that X, P coincide with the right hand sides of the continuous Neumann 
system with the Euclidean metric (|3.9I) . (|3.8p . 
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The Lax representation. As we have seen above, the solutions of the matrix 
quadratic equation (I6.4p are closely related to the Lax matrix of the continuous Neu¬ 
mann systems on Vn^r- It appears that this matrix also forms a Lax representation 
of the discrete system. 

Theorem 6.2 Up to the action of the group of reflections (13.121) . the discrete 
Neumann system (leTTD . dSJl) is equivalent to the intertwining matrix relation (dis¬ 
crete Lax pair) 

L{X)M{X, A,) = M(A, A*)L(A), (6.8) 

where 

- A)-^P A^(AI„ - A)-1 A \ 

\lr-P^{Xln-A)-^P -P'^{Xln-A)-^x)' 

^ ((A - A*)i! + r2(A*) -r(A*)) ’ 

where L{X) depends on X, P in the same way as L{X) depends on X, P and M{X, A*) 
depends on X,X in a symmetric way via (lOD . 

The proof is direct (although quite long), it uses the constraints (|3.1I) . the matrix 
identity (I3.16p and S'0(r)-momentum preservation X^P = X^P. 

Note that the 2r x 2r matrix M is a direct generalization of (|2.19p . 

It follows that, regardless to the branch of the map it preserves all the first in¬ 
tegrals of the continuous Neumann systems on Vn,r, including the non-commutative 
set given by the components of the S'0(r)-momentum 'h. Thus, when the rank of 
is maximal, generic invariant manifolds of 53^ are /-dimensional isotropic tori. 

Theorem 6.3 Let all the eigenvalues of A be distinct. Then the discrete Neumann 
system dSU) is completely integrable in the noncommutative sense with the set of 
integrals given by (j3.17p and by the components of the SO{r)-momentum mapping 

dsai). 

The proof goes along the same lines as that in [57] for case of a Liouville integrable 
correspondence, with the only difference that in our case the Lagrangian tori are 
replaced by the isotropic ones. 

Double iterations. Although the map 53^ is multi-valued, it has a kind of return 
property described by 

Theorem 6.4 Let 53“,®^ he the branches of the map corresponding to some 
arbitrary opposite partitions w_,w+. Then, for any initial {X,P) G T*Vn,r, 

<B-o®+(A,P) = (-A,-P), 

i.e., double iterations of Br with opposite partitions multiply {X,P) by —1. 
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Proof of Theorem 16.41 Expressing X from the second matrix equation in (|6.1I1 and 
using the condition X'^X = Ir, we obtain the following alternative equation for the 
matrix multiplier E, 

rur - rv - v^r - w = o, (6.9) 

\J = X^A-\X^)X, Y = X^A-^{X^)P, W = lr-P^A-\xpP 

According to Proposition 16.11 all symmetric matrix solutions of (16.91) have the form 
r = where S, T are halves of a nonspecial eigenmatrix ^ = ( E j of 


L = 




To see how the solutions of ()6.4I) and (16.9p are related, let us fix a solution T* of 
the original equation (16.4p corresponding to a certain partition w* = {rc | —w}. Let 
X*,P* be the corresponding image, U*,V*,W* be the corresponding coefficients 
in the equation (j6.9p . and {P} be the set of its solutions. Then T* G {T}. 

On the other hand, let {T} be the complete set of solution of the original matrix 
equation (16.4p but with the coefficients defined by the tilded variables X* ,P*. It is 


seen that if T 


T 


is the nonspecial eigenmatrix of 


L{XP 



corresponding to a partition {tci ,... ,Wr \ —wi ,..., —Wr}, then non¬ 

special eigenmatrix of the above matrix L corresponding to the opposite partition 
{—tci,..., — Wr I tci,..., Wr}- Indeed, 



implies 

(■w* ii) (?)—.{?)' 

for s = 1,..., r. Hence the set {T} contains the solution —T*. 

To complete the proof, consider the tilded version of the equations (16.Ih . 

P* = A^/\xpi - X*t{XP, P = -A^/\XPX* + Pt{xp. ( 6 . 10 ) 

Setting here T = —T* and assuming that P = —P, X = —X, we see that the first 
(second) equation of (|6.10l) becomes the second (first) equation in (16.11) . Hence, the 
equations (|6.10l) for X, P have the solution {—X, —P). □ 
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Real trajectories. In view of ()6.ip . for the map 55^ to be real-valued it is neces¬ 
sary that A* > oi,..., ttn, and even in this case not for any real pair {X, P) E T*V^n,r) 
the matrix equation (|6.4I) has a real solution F: as follows from the first equation 
(16.ip (the discrete Legendre transformation) and (16.21) . the real discrete momentum 
P must belong to the compact set 

X -^X (^X^A^/^{X^)X + X^A^/^{X)X'^ I X^X = X^X = l| . 

(6.11) 

Proposition 6.5 Assume A* > oi,..., an, and let the momentum Pq of the initial 
real point {Xq,Pq) E T*Vn,r belong to the set (j6.1ip . Let L(A*) be evaluated at 
{Xq,Pq). Choose the branch of the map Br corresponding to a partition 

W = {ici, ...,Wr\ -Wi, . . . , -Wr}, 

such that wi,... ,Wr are different, Wi + Wj / 0, Wi+Wj 0, i,j = 1,... ,r. Then the 
trajectory of the discrete Neumann system (j6.8p obtained from {Xq,Pq) by iterating 
the above branch is real. 

Proof. Since wi,... ,Wr are different and Wi + Wj / 0, the corresponding eigenmatrix 
'h = {ifi-'-ipr) is non-special. Then, according to Proposition 16.11 the solution 
r = T of equation (16.4p with real matrix coefficients (16.5p corresponding to 
(Xq, Pq) is symmetric. On the other hand, due to the conditions Wi -|- wj / 0, after 
setting il=U,QJ = V,2n = W, the equation (16.6p has the corresponding hermitian 
solution r obtained by the same eigenmatrix T (see Remark [7]). Since the matrix T 
is both symmetric and hermitian, it is a real symmetric matrix. □ 

Remark 8 For the case r = 1, the eigenvalues w, —w of L{X*) can be pure real or 
imaginary. The above condition gives the real trajectory for real w, as in Remark[2j 
For the case r > 1 and real A*, generically the matrix L{X*) has a set of four different 
eigenvalues invariant with respect to the inversion w ^ —w and complex conjugation 
w w. For example, consider the case r = 2 and eigenvalues w, —w,w, —w. The 
non-special eigenmatrix 'F associated to the partition w = {rci = w,W 2 = w \ 
—wi, —W 2 } dehnes the real branch of ^ 2 - 

Note that with the same initial condition {Xq, Pq) as one described in Proposition 
16.51 we can have a complex solution of (16. 4p . That is why, it is natural to consider all 
objects complexified. In the next subsection we shall give an algebraic geometrical 
description of the mapping 55^, which generalize Proposition 12.21 

6.2 Algebraic geometric description of 

Now we describe the correspondence (|6.3p under the eigenvector map 

A4 : Tred = T/S'0(r, C)/Z 2 —> Prym(5',(T) C Jac(5') 
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described in the previous section. As above, for a fixed arbitrary A* / a* choose a 
partition of eigenvalues w = {tci, ... ,Wr \ —wi, ..., —Wr} of L(A*). 

Note that the eigenvalues Wi of A(A*) and Wi of L(A*) = o(A*)L(A*) are related 
as Wi = a{X^)wi, whereas the corresponding eigenvectors can be chosen the same. 
Then let Qi = (A*, a(A*)u)i),..., = (A*, a(A*)r(;j,) be the corresponding points 

over A = A* on the regularized spectral curve S' and 

(t) 

be the corresponding non-special eigenmatrix. Let also T = T be the corre¬ 
sponding solution of the matrix quadratic equation (j6.4l) . which fixed the branch of 
the map 55^. Then T can be rewritten as 



Theorem 6.6 Under the eigenvector map A4, the above branch of Br is the trans¬ 
lation on Prym(5^,(T) by the vector given by the degree zero divisor 

T = OOi -I - -I- OOr — Ql — ■ ■ ■ — Qr ■ (6.13) 

Note that 

T -I- CrT = 2(oOi -I--I- OOr) — Ql — ■ ■ ■ — Qr — CfQi — • • • — aQr 

which is the divisor of the meromorphic function 1/(A—A*) on 5b Hence T= 0, 
and T indeed belongs to Prym(5',(T). 

We also note that the expression (I6.13p is a direct generalization of the one-point 
translation described in Proposition 12.21 for the case r = 1. In the continuous limit 
A* —)• oo we have Qi Qr '—)• ooi oOr, hence the shift vector T = A{T) 

tends to zero. 

Proof of Theorem \Q.6[ As follows from the intertwining relation (j6.8p . if '(/’(P) is an 
eigenvector of L(A), then 

is an eigenvector of L{X) with the same eigenvalue. Note that whereas ipiP) is 
normalized, the above '(/’(P) is not, so we consider normalized eignevector V’(P) = 
/~^(P)V’(P), / = (a, V'(P)) for a normalization a G 

We compare the divisors of poles of ili{P),'ijj{P) for the chosen P. First, 
note that det M(A, A*) = (A — A*)^. Hence, apart from the points of S' over A = A* 
and A = oo, M{X,X*) is non-degenerate. Therefore, apart from these points, the 
divisors of poles of '^(P) and ^PiP) coincide. 

Without loss of generality, assume that V’(P) is normalized in the same way as in 
Proposition 4.1. Then in the neighborhood of the infinite points ooi,... ,00r these 
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components have the expansions (14.111) and, near each oo^ with the local coordinate 
r = 1/\/A, the following expansion holds 


kp) 


( -r 

Vr-2l, + (r)2 +0(1)1, -Fyl V v, + 0(t) y 


( Vs + 0 (t) \ 

V'^’^Vs - TVs + 0{t)J ■ 


Therefore, in contrast to '4^{P), some components of 'ip{P) have a first order pole at 
oos- Next, observe that the eigenvectors 'tp{Qi ),..., ijjiQr) that form the eigenmatrix 
T, span the kernel of M(A*, A*), because, in view of (|6.12jl . 

M(A„ A.)T = T = 0. 


The products of M(A*,A*) with the other eigenvectors ijj{aQi), 'tp{aQr) are 
non-zero. It follows that the divisor of the above normalizing factor f{P) is 

(/) = <3l H- \- Qr + P - p - OOi - ■■■ - OOr 

for a certain effective divisor TZ. Then the divisor of poles P of V'(-P) = '4’{P)/f [P) 
equals TZ. Indeed, the zeros of 'il^{P) and f{P) at Qi + ■ ■ ■ + Qr as well as their 
poles at P -|- ooi + • • • + oo, cancel each other. Since f{P) is meromorphic on S', 
we conclude that P is equivalent to 


P + ooi + • • • + oo, — Qi — ... — 

which implies that the images of P, P in Jac(5^) differ by the translation T in 
(I6T3D. □ 


Growth of the map. Theorem 16.61 says that the shift T on Prym(5', a) does 
not depend on the step of iteration of 55,, but only on the choice of partition w = 
{wi,... ,Wr I —rci,..., —tc,}. The opposite partitions w_,w+ produce opposite 
shifts 7T,7+ such that 7T -|- 7+ = 0. Hence the 2^ branches of the map 55, and of 
the composition A4 o 55, can conditionally be divided into 2^“^ ’’forward” and 2^“^ 
’’backward” branches. Clearly, double iterations of AloiB, with opposite partitions 
give the original value of Ai{X,P). This property also holds for the map 55, itself, 
but not completely: due to Theorem 16.41 the corresponding double iterations of 55, 
give {—X,—P) and not {X,P). 

In C^, the universal covering of Prym(5^, a), all the iterations of A4 o 55, form a 
lattice H of rank < 2 ^“^. Due to the definition of the vectors A{T), for r > 2 not 
all of them are linearly independent, hence the lattice fits into a linear subspace of 
dimension less than 2'"“^. For example, let r = 3 and denote oo = ooi + 002 + 003 . 
The lattice H is generated by 4 shift vectors AI(Tl), ..., Al( 74 ) with 

71 = —Qi — Q 2 — + 00 , T 2 = —cr{Qi) — Q 2 — Qs + 00 , 

71 = -Qi - o-(Q 2 ) - Q 3 + 00 , 71 = -Qi -Q 2 - criQs) + 00 , 

and one can see that Al(71) + A{T 2 ) + -^.(Tl) + ^1(71) = 0 . 

The above implies that iterations of the ’’forward” branches of A4 o 55, give a 
linear growth of the images of the map. 
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6.3 Lagrangian description and billiards on Vn,r 

Discrete Euler-Lagrange equations. The discrete Neumann system (16.ip can 
be considered as a discrete variational problem on the Stiefel variety as well. Namely, 
let (Xfc, Pfc), /c E Z be a trajectory of the discrete Neumann system and T^, A; E Z be 
the corresponding sequence of matrix multipliers. From the hrst equation in ()6.ip 
and ()6.2p we get 

Pu = A^/\X,)Xu+i - XfeFfc, Ffc = ^ (xl^^A^/\K)Xk + XlA^/\K)Xk+i) . 

Comparing this with the second equation in (16.11) with k replaced by k — 1, one 
obtains the equations 

Xk-i + Xk-\-i = A /c E Z, (6.14) 

where Bk = Ffc-i+Ffc. On the other hand, from (16.141) and the condition X'^_^_^Xk+i = 
Ir we get the matrix equation 

Bk{XlA-\K)Xk)Bk - {Xl_^A-^l^{K)Xk) Bk - Bk {XlA-^l^{K)Xk-i) = 0, 


which determines Bk as a function of Xk-i and Xk only. In the case r = 1, when 
Bk = /3fc is a scalar, the solution of this equation has the form (I2.17p . 

Equations (j6.14p coincides with the discrete Euler-Lagrange equations on Vn,r 


dC{Xk,Xk+i)_ ^ dC{Xk-i,Xk) 


dXk 


dXk 


= XkB 


kok, 


k 


(6.15) 


of the functional S = Ylk&i, B-i^k, Xk+i) 

C:Vn,rXVn,r^R, C{Xk, Xk+l) = tT{X^A^/\X,)Xk+l) (6.16) 


first derived by Moser and Veselov in [4T]. In this sense, the expression for the 
momentum P in the system (|6.ip is actually the discrete Legendre transformation: 

r : Vn,r X Vn,r ^ T*Vn,r, Pk = - XfcFfc = A^^X,) Vfc+i - Vfc Ffc, 

where, as above, the multiplier F^ is a symmetric matrix defined from the condition 
that {Xk, Pk) belongs to T*Vn,r- In particular, the correspondence (16.3|) is symplectic 
(see [H]). 


Billiards on Stiefel manifolds We now introduce a billiard map on Stiefel va¬ 
rieties and show that its certain modihcation is a ’’square root” of the mapping 

m- 

First, consider a S'0(r)-invariant metric dsj on Vn,r defined by the Lagrangian 
Lj{X,X) = itr(V^JV). The corresponding geodesic flow is considered from a 
point of view of optimal control and discretization in [S]. There, following [4T], a 
discrete system defined by the functional S = £(Vfc, V^+i), where the La¬ 

grangian P is given by (16.160 with AAP{X^ = J have been studied. 
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Here we give another natural discrete model. Introduce the new matrix variable 
X = in which the Stiefel variety Vn^r is given by 

En^r '■ x^J~^x = I,.. 

Then the Lagrangian of the metric ds'j is given by the metric induced from the 
Euclidean metric in M"''": L{x,x) = ^tr(x^i:). In particular, for r = 1 the system 
describes the geodesic flow on the ellipsoid En^i = {(x, J~^x) = 1}. 

Consider the discrete Lagrangian 


Xfc-|-i) — |xfc+i Xfc| — ^tr((x^_|_i x/j)'^(x/j_|_i x/j)). (6.17) 

The functional S = J2f^^^C{xk, Xk+i) can be seen as a discretization of the action 
in the Maupertius-Jacobi principle. The discrete Euler-Lagrange equations read 

■j r(xfc+i Xfc) 1 Xfc_i) — J XfcCJfc, 

l^fc +1 ^k\ \^k 1 | 

where tUk is a symmetric matrix multiplier. 

Denoting 6^ = jx^+i — x^l, the system can be rewritten into the form 

Xfc+i x/j u/j J XfcWfc, (6.18) 

where 

Xfc J"^Xfc = Ir, tr(njnfc) = 1. (6.19) 

Note that the outgoing directions Vk+i satisfy the equations 

(xfc + 6kVk+i)'^J~^{xk + SkVk+i) = Ir, tr(nj+infc+i) = 1. (6.20) 

Thus, if {xk,Vk) is known, by substituting the second equation of (|6.18l) to (I6.20p . 
we get the system which determine the multipliers ((5^,0;^). 

For r = 1 this is the classical billiard within the ellipsoid For this reason 

we refer to the system ()6.18p . ()6.19p as a billiard map on the Stiefel variety En^r- 
The difficulty in studying the map, even in a totally symmetric case J = I„, is that 
the outgoing directions Vk+i belong to the variety (|6.20l) (projection of Vn,r to the 
sphere in Mn,r{^) centered at Xk), while in the case of the ellipsoidal billiard, they 
belong to the fixed unit sphere 

That is why we introduce the modified billiard map 


En..r ^ Lri, 


Er]..r ^ Vr). 


given by 

Xfc+i Xfc — Vk+i Vk — J XfcDfc, 

where now ’’outgoing directions” pk belong to the fixed Stiefel variety 

Xk J Xk = Ir, ykVk ~ ^r- 


( 6 . 21 ) 


( 6 . 22 ) 
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The symmetric rxr matrix multipliers A^, are chosen from the conditions (16.2211 . 
The mapping given by (|6.2ip . (|6.22p and (|6.18p . (|6.19p are multivalued for r > 2 
and for r = 1 they just coincide. 

One can easily observe that the modified billiard map ()6.2ip is actually a refor¬ 
mulation of the discrete Neumann system on Vn^r with ^(A*) = J written in the 
Lagrange form (I6.14p . Indeed, if we identify 


Xk = ( —Uk = ( —1)^A2A:-1, Afc = B2k+1, flfc = —B2k- 


both equations in (I6.2ip give the same equation (|6.14ll with two subsequent indices 
k. For r = 1 this relation becomes the known relation between the billiard within 
the ellipsoid En,i and the discrete Neumann system on S'^~^ described in m- 

7 Conclusion 

We have seen that the Neumann systems on Vn^r (continuous and discrete) inherit 
or naturally generalize the basic properties of the classical Neumann system on 
and, therefore, of the (odd) Jacobi-Mumford systems: the structure of the 
Lax matrices (symplectic), the spectral curve (with involution), the equations of 
motion, linearization on Abelian varieties, and, in the discrete case, the formula for 
the translation T on them. 

For this reason, we believe that the Hamiltonian systems on T*Vn,r we consider 
or, more precisely, their S'O(r)-reductions onto the Poisson manifold Q(T*Vn,r), 
represent one of the most natural matrix generalizations of the odd Jacobi-Mumford 
systems. 

On the other hand, such generalizations are quite specific since, by construction, 
the corresponding matrix residua A/) in (14.Ih have rank 1. So, it is also natural to 
consider extensions of the Neumann systems whose Lax matrices L(A) remain to be 
symplectic, but with residua of a higher rank at A = Oj. 

To reach a full generality, one can introduce the space £ of three rxr matrix 
polynomials 


U(A) = A"I^ + UiX^-^ + --- + Un, 

V(A) = PoA'‘ + HiA5-i + ... + K, 

W(A) = + WqX^ + WiX^-^ + • • • + IFn, 


where Ui,Wj are symmetric and Vq is skew-symmetric arbitrary coefficients, and 
consider a hierarchy of ffows on £ given by the Lax pairs 



In particular, N(A) can be choosen in the form similar to (|3.15p 
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which implies that Vq is a matrix first integral of the flow. One can prove that in 
the case Vq = 0 the matrix A satisfies a stationary reduction of one of the equations 
of the r X r symmetric matrix KdV hierarchy considered, in particular, in [^. 

The spectral curve S = {|IL(A) — wl 2 r\ = 0} has the involution a : (A,re) —>■ 
(A, —w) and, as above, its complete regularization S' has r infinite points. It is 
then natural to conjecture that generic complex invariant manifolds of the flows on 
£ are open subsets of non-compact extensions of the Prym varieties Prym(5',iT) C 
Jac(5'), and if Vq = 0, such manifolds are the Prym varieties themselves. 

Similarly, one can introduce a family of multi-valued maps 

Ba, : T ^ T, ([/„ Vj,Wk) ^ {ili, Vj,Wk), A* E C 


defined by the intertwining relation 


L(A) = 


/V(A) 

VW(A) 


L(A)M(A,A*) = M(A,A*)L(A), 

-v^(A)) ’ ^ ((A - A*)i, + r2 -r) ’ 


(7.1) 


where L(A) depends on the polynomials U(A), V(A), W(A) in the same way as L(A) 
depends on U(A),..., W(A), and P is an r x r symmetric matrix determined from 
the compatibility of the left- and right hand sides of (dD. This condition leads to 
the following r x r matrix quadratic equation 


ru(A*)r + rv(A*) + v^(A*)r - w(a*) = o. (7.2) 


Once one of the 2’’ solutions T of (j7.2p is fixed by using Proposition 16.11 the map 
Ba. is defined uniquely. Like in the continuous case, one finds that Vq = —V^, 
which implies preservation of the skew-symmetric matrix Vq under any branch of 
Ba.. It is expected that the branches of Ba, are given by translations on Jac(5') 
described by Theorem 16.61 

One can also introduce natural matrix analogs of even Jacobi-Mumford systems 
by constructing integrable 2r x 2r matrix generalizations of the Ragnisco systems 
given by (|2.25p and (I2.26p . 
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